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Non-Self-Adjoint Boundary Value Problems in Ordinary 
Differential Equations 


Werner Greub ' and Werner C. Rheinboldt 


(October 15, 1959 


It is shown that the theory of non-self-adjoint linear ordinary differential equations can 
be simplified and unified if, instead ot specifying linear boundary conditions in the conven- 
tional way, one merely specifies the linear subspace determined by the boundary conditions. 
If this is done, the corresponding linear space of the adjoint problem is the orthogonal 
complement of the original space with respect to the scalar product defined by the right-hand 


side of Green’s boundary formula. 


1. Introduction 


The theory of boundary-value problems for or- 
dinary differential equations of the form 
dj a, (thu (t)=f(t) (1) 


0 


is now regarded as more or less complete. The self- 
adjoint case is presented in most textbooks owing 
to its importance in applications. A complete 
treatment of the non-self-adjoint: problem is less 
common but may be found in Coddington, Levinson,’ 
and Ince. However, in these as well in the 
original publications, e.g., Bocher,’ the boundary 


as 


conditions are stated in the form 
w,(w‘a),..., a" (a) f(b), .. ., 4 (b)) 
i—1 
X , > , _ 
S* (Au (a) + Bu” (b)) =0 
(u=1, .»m) (2) 
where the @, are m linearly independent linear 
functions in the 2n dimensional number space 
RR?" (Om <2n) 
The m linear equations o,=0 (s=—!l, m) 


define a (21 —m) dimensional linear subspace 2 € [?". 
The set of solutions u(t) of the problem (1) /(2) 
remains unchanged if the functions w, are replaced 
by another system of equations @,=0 (uw=1,...,m) 
which define the same subspace 2 € P?". 


Che work was carried out while both authors were employed at the University 
of Maryland, the first in the Department of Mathematics, the second in the Insti 
tute for Fluid Dynamics and Applied Mathematics Che paper was prepared 
while both authors were affiliated with the National Bureau of Standards, the 
first working under a National Bureau of Standards contract with The American 


University. ‘Che authors’ present addresses are, respectively: Mathematics 
Institute, University of Zurich, Zurich, Switzerland; and Computing Center, 
Syracuse University, Syracuse, N.Y 
?Earl A. Coddington, Norman Levinson, Theory of ordinary differential 
equations (McGraw-Hill Book Co., Inc... New York, N.Y., 1955 
’. L. Ince, Ordinary differential equations, London (1927 
‘M Bocher, Applications ind generalizations of the conception of wdijoint 


Systems, Trans. Am. Math 14, 403 to 420 (1913). 


SO 
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This fact causes an unnecessary complication 
for the theoretical investigation of such boundary 


value problems. It results especially in a somewhat 
artificial definition of the adjoint problem. These 
difficulties do not occur if one prescribes instead of 
the boundary conditions (2) only the subspace 
Q of 7°", disregarding the manner in which Q is 
represented. If this done, the corresponding 
linear space of the adjoint problem is the orthogonal 
complement of the original space 2 with respect to a 
scalar product defined by the right-hand side of 
Green’s boundary value formula. 

It is shown in this article that the theory of the 
non-self-adjoint problems can be simplified and 
unified by a consequent use of these ideas and the 
methods of linear algebra. 


Is 


2. Green's Formula 


As usual we denote by C* (k=0,1, .) the space 
of all real-valued® functions u(t) which possess k 
continuous derivatives on a given closed interval 


la,b|. Let Liu} be the linear differential operator of 
order n(21) defined by 
L{u|=>5 a,_,(t)u” (8) (3) 
val 


C’,u(t)€ C™ and ay(t) #0 for fa, 6}. 
"~! one can compute the 


where a, _,(?) 
To every function u(t) € ¢ 
2n real numbers 

u(a).u’ (a uu" (a), ab), wu’ (b),. 2, we Pb). (4) 
We regard these 2n numbers as the components of a 
vector u of the 2n-dimensional linear number space 
7?" and call it the boundary value vector of u(t). 
The correspondence between the function u(t) 
and its boundary value vector U is a linear mapping 


Ihe restriction to real-valued functions is not essential, as can be readily 
’ 


iscertained 








u= p(u) (u(t)E Cc, 
from C"~! onto ?". 

Let 2 be a given subspace of the 7". Then the 
boundary value problem in consideration assumes 
the form 


Liu|=/(t), 


p(we Q, (wx) 


with given /(f)¢ ©°. The corresponding homoge- 
neous problem is 


Llu|=—0, 


p(uye Q. (Wz) 


For any two functions u(f), e(f) © C", the integral 


*) 


(Ll ul], v) L\u\ edt (6) 


aa 


is defined and by repeated partial integrations can 
be written as 


(Li{u], v)=(u, L*[e]) +[u,r]i, (7) 
where 
L*(v]=>5(—1)"(a,-») “ (8) 
1 i 
n—1 n 
[u,v]2=S>) >> (—1)4-"" "ua (a,_ we '|2. (9) 
v=0 p=v+l 


Equation (7) is the well-known Green’s formula. 

The linear differential operator L* is called the 
adjoint operator of L. The expression [w,v|) is a 
bilinear form in the two boundary value vectors 
p(u), p(v). For arbitrary vectors 


x Se y aie We ona t cag ie 
of the 722", this bilinear form can be written in the 
form 
n—l a~e—i a a x 
gxy=>> DS Dd (-1)(") 
v=0 yp=0 =p re 


[a5 -1(b) Fen —a 25”, -1 (a) EO n(?’]. (10) 


With this notation we get Green’s formula, 


(Liu), D) (u, L*[v}) tyl pl), p(r)) (11) 
It is obvious that the correspondence L —> L* is 
linear. Furthermore, the adjoint operator to L* is 
again L. ‘To show this, apply Green’s formula to 
any ue C™ with p(w) =pl(v)=0: 
(L{u], vp) -(U, L*[v)) - 0 
(L*[v], uw) —(v, L**[u])=0, 
which results in 


(L 





uj]— L**[u], v)=0 
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(for all wee C” 
to 


with p(w p(v)=0) and thus lead 


L\uj—L**{u) 
for all functions u(f)¢ ©". In the last step one has 
to use the same technique used in calculus of varia. 
tions for the proof of Euler’s condition. 

The differential operator L is said to be self-ad join 
if L*=L. On comparing the forms (3) and (8) gf! 
L and L* one sees immediately that » must be eye 
if the operator is self-adjoint. 

In the case of a self-adjoint operator L the left 
hand side of (11) is skew-symmetric with respeet 
uw and v, and this implies that the bilinear forp 
¢g(U,V) is also skew-symmetric. 

We will now show that ¢(X,y) is not degenerate. 
To prove this it is sufficient to verify the following 
statement: If for a fixed vector ye 2°" we have 


¢(X,y)=0 for every xG 2°", then it follows that 
y=0. 
We prove this first for a fixed vector y of the form 
(a) a ' 
Vo No > Wn—1» 0), aca ee 0 (12a 
or 
y,— {0, aes (12b 


taking s correspondingly as x, or X,, analogously de- 


fined. In this case the expression g(Xo, Ya) for a=a 
or a= is a bilinear form in the two n#-dimensiona 
vectors £0 &e)) and n> ey 


Therefore, if ¢(Xa, Ya) —0 holds for every x, and 4 
fixed y., then in expression (10) the coefficient of 


every &2 must vanish separately. Thus 
n—v—1 n—v—!1 
‘ ‘ 5 
24 2u (~1) CG) e%.2,-: @ 99=0 
B= A= 
(y=0 ,n—!l 


For y=n—1 this has the form ay(a 
a,(a) #0 it follows that 


nN (). 
For n=n—2 we get now 
a,(a@) ne ao (a No Ay(a@) ny Ayla) n'y 0) 
and 
nN) 0. 


Continuing in this way, we ultimately show that 


ny = (u=0, . . ., n—1) 
and thus yz=0. 

Now the proof is readily completed. 
tion (10) it is seen that for arbitrary x=x,-+ 


Y=YatY>, 


From defini- 
x, and 


o(X, Y)=¢(X», Yn) — (Xa, Ya). 
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For a fixed y this expression shall be Zero for every 
x C2" and in particular for arbitrary vectors x=x, 
or X=Xp. Hence the two terms g(X,, y,) and g(X», Y») 
have to vanish separately for every x, and X,, re- 
spectively. In accordance with the first part of the 
proof, this Is only possible for y.—0 and y,—0, Le., 


0. 
. 3. The Adjoint Problem 


The bilinear form ¢ was shown to be nondegener- 
ate. Accordingly, it can be employed to define 
“orthogonality”? of vectors x, ye /?". We will say 
x and y are g-orthogonal if g(x, y)=0. It should be 
observed that this g-orthogonality is not necessarily 
symmetrical. 

“One can speak, in particular, of the g-orthogonal 
space to 2 that is the space Q+ of all vectors yC R°’ 
which are g-orthogonal to all vectors x€Q. 

For the dimension of Q+ evidently holds 


dim Q 2n—dim Q=m. (13) 
Equation (13) implies 
dim (Q2 dim @, 
and henee, together with (2Q+)46€Q, that 
(G3+)- Q. (14) 


It is worth while pointing out that the relation 
Qn2+—0 is not, in general, valid; there may exist a 
nonzero vector x©Q2 which its g-orthogonal to all 
vectors of Q. 


With this terminology we define the adjoint 
boundary value problem of (ay): 
(wx) L*\uj=g(t), p(u) EQ 
with given g(f) “hes 


From L**=—Z and (14), it follows that the adjoint 
problem of (#}) is again (my). 
In particular, (ay) 


L*=L, Q—04 


(] 5) 
4. Homogeneous Problems 
the 


Consider homogeneous problems (a,) and 


th). It is known from the theory of homogeneous 
linear differential equations that all solutions u € CC" 
of the equation 

L\u|]=0 (16) 


form a n-dimensional function space A. The solu- 
tions cr CO" of L*|r]=0 form another n-dimensional 
function space A*. These two spaces A and K* 
are mapped by the operator p into two subspaces 
K=p(A) and K*= p(K*) of PR". 


* For self-adjoint problems Green's formula states 





Liu], t (u, L{r] for Q 


plu), pit € 


Usu lly one uses this relation to define self-adjointness of (wy The two defini- 
tions are equivalent as can be shown in a way similar to that used for the proof of 
L*=L, 


will be termed self-adjoint ° if 


In these subspaces K and K* the operator p has a 
unique inverse p™'; hence, A=p-'(K) and ix’ 
p ‘(K*). This is immediately clear because p(u)=0 
implies that all initial values w(a) (v=0,1, . . ., 
n—1) are zero; hence u(t)=0 owing to the unique- 
ness theorem for the initial value problem of eq (16). 

The operator Z and the bilinear form g(x,y) deter- 
mine the ¢-orthogonal space K= p(A)+ of K. 


THeorem 1: The subspaces p(KX*) and K+=p(Kk)+ 
of RR?" are identical, 1.¢., 

K+ = p(K*). (17) 

Proor: For any two functions u(t)€ AK, v(t) C A* 


Green's formula gives ¢(p(u), p(v))=0, which implies 


p(K*) € p(K)+=K4 (18) 


Since the mapping p is regular in AK and K*, it follows 
on the other hand that 


dim K 


dim p(y) n 
and 
dim p(A*)=dim K*=n. 
Thus, we find 
dim K+—2n—dim K=n, 


and the inclusion (18) must constitute equality. 

To this point the boundary conditions [in (ay) and 
(77;) have not been used. ‘To consider them now we 
introduce the space A of all admissible functions of 
(7), 1.e., the space of all functions u(t)€ C" with 
p(uje2. Thus, the solutions of (7) are, exactly the 
functions 


u(t) C ANRK, 


(19a) 
For these functions the boundary value vectors are 


p(uy CON p(K),. (19b) 

We handle (7%) in the same way. Let A* be the 

space of all 7(t)€ C" for which p(w) € 2. Then, the 

solutions of (77) are exactly the functions 

v(t) EC A*NA*. (20a) 

In view of relation (17) the boundary value vectors 
of these functions are 


e(v) € 24NK+A= 04Np(kK)-+. (20b) 


Thus, the number of linearly independent solutions 
of (7) or (7%) is given by 
k=dim (ANK) 


dim (QN p(4)) (21) 


and 


| 
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k*=dim (AN A*)=dim (Q4Np(K)+), — (22) 


respect ively. 








THEOREM 2: 


The numbers k and k* of linearly | 


independent solutions of (aq) and (xj) are connected 


by the relation 


k*—m—n-+k. (23) 


To prove this we use the following well-known 
lemma from linear algebra: 


Let A be a linear Space and Aj, A, two subs pace 8, 


then 


dim A,+dim A,=dim(.A,0.A.)+dim (4,+ A, 
(24) 
where A, + Ag is the linear closure of A, and As. 

This relation is employed in a slightly different 
form. We introduce the orthogonal subspaces 
AEA and AFEA of A, and Ay with respect to 
a given (nondegenerate) bilinear form. Whence, 
we have 


dim Ay=n—dim A, 
and 
dim (A;NAs)=dim (A,+ A,)+=n—dim (A,+ A?). 
Thus (24) becomes 
dim A,—dim (.A,Q44.) =dim Ay;—dim (AtNAS). 
(25 
In the present case, 
A=", A,=Q, 1.—p(K 
therefore 
dim A,=2n—m, dim Ay;=m, dim uls=n;: 
hence, 
n—dim (QN p(Ay)) =m—dim (24M p(Wy)+), 
which completes the proof. 
5. The Nonhomogeneous Problem 
Consider the nonhomogeneous problem (my) to- 


gether with its adjoint homogeneous problem (7/ 

We shall prove the “alternative” theorem of the 
theory of systems of linear equations. 

THEOREM 3 (Alternative theorem): (ry) has a solu- 
tion if and only if (fv) 
(ri). This can be expressed in the form ‘ 

Li A|=(A*NK*) (26) 
denotes orthogonality in the 
with respect to the sealar product 


Here the symbol 
function space C° 
(f,q). 


? L[A] is the subspace of C° containing all functions L[u] with w€ A. 


0 for all solutions v(t) of 
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Proor: (1) Assume (my) is solvable and u(t) is q 
solution. Then, for any solution r(t) of (x), it fol. 
lows from Green’s formula that 
(u,L*\r}) 


(2) ), 


(fv) =(L[u],v) t+ (p(u),p 
Since L*{r|=0, then (u,L* le} 
and p(v) € 2+ it follows that ¢g(p(w), p(r)) 
it is seen that (f,v)=0. 

(2) Assume that (f,7)=0 for all solutions of (x), 
We shall show that (r,) has a solution. 

From the existence theorem for ordinary differentia] 
equations it is known that the differential equation 
Llu| f must have at least one solution u,(f) 
For this function u(t) and any solution v(t) of (x3) 
we find from Green’s formula that 


0; and from p(u) €Q 
0. Thus 


0) (Liu), ) ( Qo. L*\r 


4 pir 


(f, i) ol ply), p(v) 


yl p| ui 
Thus p(w) is g-orthogonal to all p(r 21179 p(K)+ 
and 
QO+ p(y 


p(y) e(Q p(t 


where the sign again denotes the linear closure 
of 2 and p(y 

This inclusion states the existence of two vectors 
u,€2, we pl(A) such that p(m)—u,+u,.. From the 
fact that p has a unique inverse p | in K it follows 





that there is exactly one function u,(f p '‘(w)CK 
for which L{w.|=0 and plw)=—u.. We form a new 
function u(t) %(ti)—u(t This function u(t) is a 


for obviously Lluj|—/f and plu 


(>) 


solution of (my, 

p(y) p(s) u, 
Taking the orthogonal complement on both sides 

of (26) we find 

A*NK*. (27 


LiA (A*NK* 


Thus we get the following 


Corotiary: Lf for a given function v(t) &¢ 


(v(t) Lhu (0) u(t) eA, 





for every 


then L*¥{r]|=0, pl(v) © 2 and in particular r(t) EC". 


6. Green's Function 


The solution of the nonhomogeneous problem 
(ry) provided that it solvable determined 
up toa solution of the homogeneous problem (Tq). 
In order to obtain an explicit formula for the solution 
of (my), it necessary assume an additional 
condition which uniquely determines the solution 
of (ry). Furthermore, it is desirable to modify (ty) 
in such a way that it is always solvable. 

With these two aims in mind we introduce an 
orthonormalized basis in each of the spaces 
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and A*N K*. 
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Let 
,, Ue, u, and Dis Oa « -— 2 
be these two bases. With the functions ¢,(t) we 
define the new boundary value problem 
Liu} f(j—>> (f(t) w(t) )er(t), p(w) € Q). (29) 


Owing to the orthonormality of the v(t), the right- 
hand side of the diffe rential equation is ‘orthogonal 
to all v(t) (v=1, . _ k), and thus to all ‘solu- 
tions of (77). There fore (29) is solvable for every 
function f(t) C°. To find an additional condition 
which uniquely determines the solution of problem 
(29). we recall the following fact: If w(?) is a solution 
of (29), then 


k 
u(t) +>) y,u,(t) 
l 


Me 
is also a solution for any real coefficients y,. We 
can determine the y, in such a way that 
(p==1, =. 


( u(t) +>) yuuy(t), u(t) ) 0) 
w= 


Owing to the orthonormality of the u,, these rela- 
tions are satisfied for 
Vu (u(t),u,(t)) (n=1, op EDs 
This shows that there is always a solution of (29), 
which is orthogonal to all uw, and thus to the whole 
space ANA; obviously, this solution is uniquely 
determined. 
Consider instead of (xy) the following problem: 


k* 
L\u\j=f(t)— S35 (f.x,) v(t) 

u=l 
p(uyEQ (wy) 
(u(t),u,(t))=0 (u=1, , &). 


We know that this problem (zy) has exactly one 
solution u(t)€ (" for every function f(t)E€C°. The 
correspondence between u(t) and f(t) is denoted by 
an operator G, L.e., 

u(t) =G[f(t)). (30) 
This operator is obviously linear. It will be shown 
that G@ is an integral operator of the form 


Dd 


q(t,r) f(r) dr, 


a 


u(t)=G|f (t)] | (31) 


where g(t, 7) is the so-called Green’s function in the 
generalized sense. 

To prove this, we consider functions g(t, 7) having 
the following properties: 

(P 1) g(t, r) is a continuous function of t, 7 for 
asSt.r<b; 
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(P 2) For every arbitrary but fixed r=79,a S77) Sb, 


g(t,7>) possesses ” continuous derivatives for 
astSry and 7 Sth such that 
0" (To T 0, To) Oo": ( 0 (), o) 

: one ™’=0 (v=0,1, —2) 

ot’ Of, 
oO” la (To +-Q, To) oO” 'g(t>—0, To) I - 
or ot"! “Ao (To)’ 

(P 3)° For every arbitrary but fixed r=r) in 

asrmsb, 
k* : 
Li g(t, t.}=— DS oy (t) vp (0) 
pel 
fora stsb but t¥ 7: 
(P 4)° For every arbitrary but fixed r=7r)ina<19)<b, 
0: (4 (t,po)) €2: 

(P 5) For all functions u,(t) (w=1,.. ., k) of 

the basis (28) 
*> 
| q(t,7) Uy (t) dt=0 
Ja 

identically in arb. 

It shall be proved that there exists exactly 


one function g(t,r) having these properties which 
satisfies (31). 


me prensa bo (P 1) “inane (Ps 5), it is “uniquely 


determined. 


Proor: Let g(t, 7) and g(t, tr) be two functions 
which both possess properties (P 1) through (P 5). 
For an arbitrary but fixed r=7) in a rb we define 
the function 


go(t) = gilt, To) — galt, To). 


Then, according to (P 3) it is seen that 


L\go(t)|=0 for astsb, t#rm. (32) 
In view of (P 2), go(t) € C"'; hence it follows from(32) 
that the jump in the (n—1)- th derivative of go(t) at 
t=r) must be zero. Thus, go(t)€ C"~' and in view of 
(32), go(t)€ K. From (P 4) it follows that go(t) € A; 
therefore, g(t)€ ANA. On the other hand, it can 
be concluded from (P 5) that go(t) € (ANA)"*. This is 
possible only if go(t)=0 for a <tx<b. Bee ‘ause T>) Was 
arbitrarily chosen in the open interval a< r< 6, this 
is oquavarent, | to the statement g,(t, 7) =golt, r) for 
ast<b, a<r<b, and in view of (P 1) for the whole 


ast,rsb. 


* The subscripts ¢ in L:{g] and p:(g) denote that the differential operator Z and 
the mapping p operate on g(t, r) regarded as a function of ¢. 








To state the existence of g(/, 7) we first prove the 
following: 

Lemma 1: There exists at least one function h(t, r) 
which POSSESSES prope rtie Ss (P 1). (P 2). (P 3). and 
(P 4). 

Proor: For fixed arbitrary 7 in a— r~ 4, 
the ordinary differential equation 


consider 


k* 
L, (g(t, r)|=—>S v(Oe,(7) (33) 
p= l 


and the two sets of initial values 


O’g(t, rT) 


0 (yv—0, 1 n—1) (34a) 
or’ — 
and 
oO” !a(t.r) | 
J , (34b) 
or” sme Go(T) 


Both initial value problems (33) /(34a) and (33)/(34b) 
possess unique solutions h,(f, 7) and ho(t, 7), respee- 
tively. These solutions A; and fy are continuous 
functions of ¢, + in the whole square, ast.rsh. 
Furthermore, they possess » continuous derivatives 
with respect to ¢ in this square. Aecordingly, the 
function 


h,(t,r) for estsr 

ha(t,r) 

h(t,r) for r<tsh 
obviously possesses properties (P 1), (P 2), and 
(P 3). 

We define the coefficient-functions 
Cy(T) (ha(t, +), u,(t)) (oerd. . sng RD 


which in view of the continuity properties of h(t, 7) 
are continuous functions of 7. Hence, the function 


Ss ¢,(7)u,(t) 
p=l 


h(t.r)=he(t,r) 


evidently satisfies properties (P 1), (P 2), and (P 3 
Furthermore (P 5) obtains. since 


(As(t, T), Uy(t))—e,(7)=—0 


(ul sia, 


(A(t, r), u,(t)) 


owing to the orthonormality of the u,(¢). 


LemMaA 2: There exists a_ linear transformation 
which maps every vector XCQ+ pK) into one vector 
y € p(K) such that 

(p(y), u,(t)) =0 (g=rl, .. . . &). (35) 

Proor: Let x€Q+ p(K) be given. Then there 

exist at least two vectors X,,X, for which 
x=—x,+X, x, 62 x, € p(K). 


x, and xX, are determined to within a vector belonging 
to the intersection QUp(/Ay). With coefficients 


Cu=(p '(Xe), u,(t)) To) ees | 


the following new vectors are defined: 


k j 


yi=%X1i74 Cue | Uy), Yo=X SSe,p(u,). 


; p=l u=l 
since 
dSc,e(U,) ON pk), 
“ l 
then 


X=y,+y, y€2 y> p(k), (36 


and, furthermore, in view of the 
the Uy, 


orthonormality of 


k 
Sc (p ‘o(u,), uy) 
I 


(p~'(X.), uy) —c,=—0 (u=1,.., &). (37 


(p '(¥2), Uy) =(p 


The two vectors y, and y, are uniquely determined 
by conditions (36) and (37) This is rather obvious: 
for any other decomposition x=—y;+yYy, of xX which 
satisfies (36) and (37) it follows that p~'(y.— ys) © ANK. 
However, the relations 


, 


(o-'(yo—y3), u,)=0 (u=1,...,&) 
obtain only if y,=y, which then implies y,=y}. 

Hence, it has been shown that we can 
with every XC 2+ p(y) exactly one vector y. 6 p(K 
which fulfills (36) and (37 This was the statement 
of the lemma. 

It is necessary at this point to obtain an extension 
of Green’s formula for functions having discontinuity 
properties similar to (P 2). Let the funetion y(t 
possess the following properties: 

1) y(t) continuous derivatives for 


i possesses NM 
astsSr, rSts<b, where rc is a certain fixed point. 


~ 
to 


where ¢ is a certain fixed constant. 
For any arbitrary chosen r(f) © C", it 
shown by repeated partial integration that 


can be 


(38 


CL 


DIT). 





y\.v) (y,L*[r]) pir C Ag T 


Y\ ply), 


This is the generalized Green's formula. 


There exists a function g(t,r) whieh 
1) through (P 4). 


Proor: Let A(t,7) be the function which has prop- 
erties (P 1), (P 2), P 3), and (P 5). For every 
arbitrary but fixed 7 in ac 7b, insert A(t,z) and 
the basie funetions »,(f) (v=1, .. ., &) in the 
generalized Green's formula: 


THEOREM 5: 
POSSESSES the five prope rties (P 


(L, lA}, vy) (h, L*\r,| 


t o(p(h).plr, ) ry(T 


go(p,(h), p\ry) VyCT (T i. o-<% k*). 


associate | 


On 
(P 


an 


Wi 


{10 


~ 


be 
{10 
a . 
col 
pr\ 
ol. 


It 
ful 
de! 


TI 
(P 


Fi 


ging 


(36 


ity of 


(37 


ined 
V10us: 
Which 


ANK. 


Yi. 
erate | 
p(k) 


‘ment 


nsion 
nuity 
1 y(t) 


s for 
it. 


— 


n be 


rh ich 


rop- 
very 
and 
the 


k*). 


function 


On the other hand, for ast<6, t#r, according to 


(P 3), 

~~ 
— 
p=l 


(y(t), v(t) )v,(7) v,(7), 


(L,[h], v) 


and in view of the orthonormality of the v,(f) 


o(p,(h),o(v,)) 0) fou: f. , &*). 
This is equivalent to 
p(h) € (24Up(K)+)4+=2+ p(k). 


We now consider the (continuous) linear transforma- 
tion from 2+ p(/x) into p(X) which exists by lemma 
9 Let h€ p(y) with 

(p'(h), u,)=0 (a=1, _k) (39) 
be the image vector of p,(h) under this transforma- 
tion. Then h continuous function of + for 
axrsb. This is obvious because p,(h) depends 
continuously on 7, and the correspondence between 
o(h) and his continuous. Let x,,..., X, be a basis 


of p( AX) and 


IS a 


" 
h=>> d,(7)x,. 


w=l 


It is clear that the coefficients d,(r) are continuous 
functions of T for as rah. Hence. the (uniquely 
determined) function 


" 
be, 
—_—- 
“ 1 


h(t. 7) p ‘(h) d,(r), p '(x,) K 


iscontinuous in both variables ¢ and 7 for aSt,rsh. 
We define the function 


q(t,r) =h(t,r) —h(t,r). 
Then g(t,7) obviously possesses properties (P 1 
(P 2), and (P38). (P 4) follows from 
p:(g)=—p,(h)-he€ Q. 


Finally, as a consequence of (39), 


(q(t,r) u(t) (h(t,r) u(t) ) (h(t,r) u,(t)) 0) 


since A(t,r) fulfills (P 5). 
(P1) 


for s=1,.... k, 
Thus 
(P 5). 
Now we are in a position to prove the following 
final theorem: 


q(t,r) properties through 


This proves the theorem. 


pe SSESSES 


Let g(t,r) be the uniquely determined 
the (P1) through (P5); 


THEOREM 6: 
with prope rties 


then 


u(t)=Gl[f(t)] g(t, rt) f(r) dr, 


Ja 
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for every given f(7) © C°, constitutes the (uniquely de- 


termined) solution of (ty). 

Proor: Let v(t)€ A* be an arbitrary function. 
For any fixed r in a<_r<b we apply the generalized 
Green’s formula (38) to g(t,7) and v(t). With regard 
to (P 4) and p(v) € 2+ it is seen that 


(Lilg], ce) =(g, L*[v])—v(7); 


thus, in view of (P3), 
k* 


2 (0s) v, (7) 


(L,l|g], v) (g,L*{[v]) —v(r) 


for astsSb, a<r<b, t=r. We multiply this 
equation by a continuous function f(r) and integrate 
over r. Because g is a continuous function of ¢ and 
r, it follows that 


k* 
‘ > . - e 
2a (Pte) (of) 


(Gf), L*[v]) — Jr). (40) 


For f(t) we consider the boundary value problem (29) 
which is solvable no matter how f(t) was chosen. 
Let u,(t) be a solution. For the given function 
v(t) A* we obtain 

- 
(fv) — DSF e.)(e%,2), 

v=l 


(Li u,},v) 


and thus together with (40) 


(Llu), 0) = (GLfl, L*tel). 


Applying Green’s formula (11) to the left-hand 
side we conclude that 


(L* lr], uy) = (GLP), L*lr])—e(e(u), e(v)), 


and because g(p(u,), p (v))=9, then 


(G[fJ—m, L*{r))=0. 


This relation is valid for every function v(t) € A*. 
Hence, it follows from the corollary to the alterna- 


tive theorem that G{[f]—u,€C"; thus, u(t)=6[f] 
Cc”. Furthermore, it follows that 
L\G|f|—u,]=09, p(G| f]—u,) €Q, 
which is equivalent to 
LGN = Llu) =f—S.0) 


o(G{f]) € 2. 


Finally, it follows from (P5) and (P1) that 


(u(t), u,(t))=0 (g= 1, os &). 


In applying theorem 6 to the adjoint problem 
(#y*) we confirm that for any f(t) € C° the uniquely 





determined solution v(t) of (#x*) is given by | With another continuous function f(t) we find frop 
Green’s formula (11), 





*b 
v(t) =G* [fi] | g* (tr) f'(r)dr. (41) . a 
ai (G|f\f,) (GP LAIG*| f, )- SS(fi.u \(G|f\ uy) 
v= 
We will show that the integral operator G* is the a cae oa iat i cet he 
adjoint operator of G, Le. 7 : (LIG(F\,G@* fi) + e(o(Glf)),0(@*| fi) 
(GUL A=, GLAD (42) Do (fim (Gf | ,ue) ' 


4 te Co . a or er 
for all f,,/ ° (f,G*{f,]) . fr.) (v G*| f;\) 


Let v(t) =G*[f\], ie., —' ws u 
k lle -_ : : 
L*(@*|fill=h Dhue) Ur tel p(G| sf |), e(G@*| fi))) 4 Do SiMe) (GLP) ,u ; 
v=l1 l 
*[f,]) €Q4; (G*{filed=0 | f, wiag ee ¥: — — , » oe : 
AG" fil) € Pl Sad ee “ Since G[f,] and G*|f,] are solutions of (a#y) and (#,*), 
respectively, the last three quantities must vanish 
which leads to (42). 
From (42) follows for the corresponding Green's 
functions, 
Wasuinoton, D.C. (Paper 64B2-24) | y*(t,7) =g(1,t). 
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Criteria for the Existence and Equioscillation of Best 


Tchebycheff Approximations 


John R. Rice ' 


(January 


a 
}},u Let F(a,r) = Za,e,(r) where 

The following theorem answers the 

of Tehebycheff approximations: THreorem. Let 
(ty*). f(r) be an arbitrary function continuous on [0,1]. 
“anish mation, (B) a necessary and sufficient condition that 


f(x) is that max|F(a*,xr) —f(a) 


alternates at least » times in [0,1], 


9, 1960) 


¢i(4) are continuous in [0,1! and the a; are real parameters. 
principal questions of a general nature 


in the theory 
form a Techebycheff set and let 
Then (A) f() possesses a best approxi- 
F(a*,«) be a best approximation to 
(C) the best approxima- 


bei()} 


— tion to f(%) is upique. A. Haar [Math. Ann. 78, 43-56 (1928)] posed and answered the 
een s following question: What conditions on F are necessaiy and sufficient for theorem C to be 
valid? The condition he found is that {¢,(7)} must form a Tchebycheff set. This paper 
poses and answers the following three questions: (1) What conditions on F are necessary 
and sufficient for theorem A to be valid? (2) What conditions op F are necessary and 
' sufficient for theorem B to be valid? (3) What conditions on F are necessary and sufficient 
for both theorems A and B to be valid? 
This paper does not tacitly assume that the a; may assume all values. 


|. Let 7 be a set of n functions g(r), 
y, (4) continuous on [0,1] and let F(a,r) 
be called a 7-polynomial. F(a,r) is said to be a non- 
trivial T-polynomial if Sft.yja,;\>0. It is assumed 
that there is no nontrivial 7-polynomial which is 
identically zero. The parameter set P of F is the 
domain of the a,. 7 is said to be a Tehebycheff set if 
every nontrivial 7-polynomial has at most n— 1 zeros. 
max! F(a,r)—f(r)| is said to alternate n times if there 
are n+ 1 points 0<.;- Snail such that 


I | 


max F(a,r)- f(x) 1 


go\r), 
p> - ja (0) 


- 


Ip: 


F(a,r;)—f(2r;) [F(a.r,.,)—f(a, 


All maximums and minimums are taken over s€{0,1] 
unless otherwise stated, and all summations are from 
! to n unless otherwise stated. Braces, denote 
sets and sequences, and (4 read “the 
set of x such that 

The following theorem answers the principal 
questions of a general nature in the theory of Tcheby- 
cheff approximations. 

TueoreM: Let T be a set of n functions and let f(x) 
be an yo a function continuous on [0,1). Then 

A. f(r) possesses a best approximation, 

B. a necessary and — tent condition ey F(a* x) 
be a bestapproximation to f(r) is that maz| F(a* xr) —f(2)| 
alternates at least n times. 

(. the best approximation to f(x) is unique. 

A. Haar? has posed and answered the following 
question: What conditions on F' are necessary and 
sufficient in order that theorem C be valid for every 
continuous f(r)? The condition he found is that 7 
must be a Tchebycheff set. 


ee 
' Part of this work was done while the author was an NRC 
e8e arcn Associate 
7A. Haar, Die Minkowskische Geometrie 
je n, Math. Ann. 78, 43-56 (1928). 
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Three similar questions are posed in this paper. 
They are: 

(1) What conditions on F are necessary and 
sufficient for theorem A to be valid for every con- 
tinuous f(r)? 

(2) What conditions on F are necessary and 
sufficient for both theorems A and B to be valid for 
every continuous f(z)? 

(3) What conditions on F are necessary and 
sufficient for theorem B to be valid for every con- 
tinuous f(x)? 


These questions are answered by theorems 1, 2, and 
3, respectively. 

In this study it has not been assumed that ?=£,, 
and indeed the topological nature of P plays a 


central role in the analysis. 


The following definition will be useful: 

Siemamaneae F has property Z if a¥a* implies 
that F(a,r)—F(a*,x) has at most n—1 zeros in {0,1}. 

Note that this does not imply that 7’ is a Tcheby- 
cheff set for the parameters are restricted to P. 

The first theorem gives an answer to the first 
question posed. 

THEOREM I: 
closed. 


PROOF: 


Theorem A is valid if and only if P is 


Let d=inf max |F(a,r)—f(z)|, then it can 


acP 
be shown * that there is an M< @ te that \a,|>M 
implies max) F'(a,r)—f(x)| >d- f(x) has a best 
approximation if there is an a* ep such that 


max | F'(a*,2)—f(r)|=inf max | F'(a,r)—f (2)}. 


Let 
Py 


’N. I. Achieser, Theory of approximation, p. 
New York, N.Y., 1956). 


fala€P,>3\a,|<M}. 


10 (Ungar Publishing Co.» 





Clearly a may be restricted to values from Py. If P For the purposes of the next lemma an arial gy 
is closed, then inf max! F(a,r)—/f(x)! is attained and | will be defined. A set in /, is said to be axial jf 


a€Pu every neighborhood of each point @ of the set cop. — se! 
J (x) possesses a best approximation. tains a point of the set on each half-axis (positive ay th 
Assume that P is not closed and let a’ be a bound- | yjegative) of the /, coordinate svstem with origin CO 
ary point of P not contained in ?. Consider the | translated to a. “a TT 
function Sajg,(r). If theorem <A is valid, then it Since the validity of theorem B implies that pf sit 
has a best approximation, say F(a*,r). Now it is | has property Z, the parameter space P may y If 
clear that identified with the function values of F at n fixed? th 
inf max |! F(a,r)—>3 aig; (4) =0, distinct points in {0,1| Let Om rs r, a ts 
acP r,—1 be n such points ‘and let P’” denote the set pl: 
hence identified with P. an 
» (a; —a¥) ¢, (4) =0 for ré [0,1] LeMMA 2: If theorem B is valid and P’ is closed. (| 
eae yer ; : _— then P’ is an arial set. es 
This possibility is precluded by one of the initial Proor: Let €>0 and a point b*€P” be given 
ASSUMPLIONS ; SO a contradiction has been reached. F(b*.x) shall denote F(a*.r) where a* € P is identified) Lé 
Phe following example shows that the assumption | ith ’,* — Let | be 
that no linear combination of the g,(2) be identically Ties m= 
zero is necessary for the validity of theorem 1. Let d=min [2,,02—21, - - +5 4n—Lna1,1 —2p). \, 
ei(r4)=1, g(r)—1, and let P be the subset of F, | of 
consisting of ((7,0)\r4#1} and (0,1). Every con- | A sequence of functions | f,,(2)\m—2.3, -}, each is 
tinuous function has a unique best approximation, | continuous on [0,1], will be defined as follows: on 


but 7? is not closed. , ' .. 
Ain (x) = F(b* 2) + (—1)’e/3) for rE [O,r,—d/m), § ist 


3. In order to answer the second question, three 





. ° ry. . e - > | > | ~ » a ¢ 
lemmas will be established. The first is rE [ajtd/m,rj,,—d/m], j=1,2, . . .,k-2;) be 
Lemma 1: Jf theorem B is valid, then F has — _ 
property Z. | Im(4) = F(b* a) + (— 1) 'e/3 for if 
) me _ 7 : PYia® a ae a oF 
Proor: Assume F(a,r)—F(a*,r) has n zeros. Let | x € [x,_,+d/m,r,,,—d/al: . 
M(r)=3|F (a,r) + F(a* - 
- , “he 7. ’ a 
Am (2) = F(b* a) + (—1)?*"e/3 for rE [r,+d/m,1], 4 
d=} max F(a,r)—F(a* ap ly 
| r€[rj,+d/m,rj,,—d/m|, jok+l, .. 0 
and let O<2,<72<, .. ..<2na,;<1 be n pomts 
y re I” = * » ‘ _ . . : . . ° 
where F(a,r) Fa vv 0, along with One pot Fir) IS defined in the re maiming intervals so that 
, re . — he , * » fe . ° ’ » : : 
where F(a,r)—F(a*,1) 40, say F(a*,r,)>F(a.rz) | | f,,(2)—F(b*,2)\ 3 and f, (2) is continuous. in 
for concreteness. Set (0,1). 
= Since P”’ is closed, ? is closed and every function 
6=% MIN (4), F2—-T}, ~« ~ +, Engi —Ln, 1—Fn41), has a best approximation by lemma 1 and theorem 1. 
uy Let F(b,,0©) be the best approximation to f,(z 
with yin or ] — Fn l omitted if they are Zero. Now it is clear that sen | F (h,, a) ) I (b* x, | (— ] A ! 


J 





A function f(x) continuous on [0,1 
as follows: 


will be defined | and that F(b,,.2) —F(b* 7) has a zero in 


f(r)=M(r) in the intervals” [z,+6, x,,,—8}, [4j—d/m, 2, +d/m)\,j—1,2,..,k#-1A +1, he 
geil 2... 2} Since ?” is closed (and clearly the 6,, are bounded), 
F : es b,, contains a convergent subsequence with a limit 
I2)=M(z) in [0,2,—6] if >0; in P’. Let by be this limit. Now F(b*,2)- 
f(z)=M(x) in [aea,+6,1) if aya,<1; F(by.7;)=0, 7 - eer ks 1e+1,..., ne By 
; ae a theorem B, max | F(6,,.1)—f,,(4) alternates n times; 
in the remaining intervals f(r) satisfies | hence 
! , . 
f(x;) = F(a*,x;) + (—1)?*"(3d/2) max Fb) — F(b* 0) >€/3, 
and , ’ 
since max |F(b*,7)—F (2) alternates only n—l 
|F(a,r) —f(r)|\< 3d/2, |F(a*,2)—f(r)|< 3d/2 for rx, times. Therefore, 


T , , max |F(b)1)— F(b* 7) > €/3 and 6, #b* 
Now max |F(a*,rc)—f(r)) alternates exactly 7 — 


times and max | F(a,r)—f(r)| alternates exactly n—2 fin 
times. By theorem B, F(a*,r) is a best approxi- 
mation to f(r), and since 


v actual construction it has been shown that 

there isa by€ P’ such that dbp is on the kth coordinate 

axis (positive or negative according to whether t 

max |F(a*.r7)—f(r)}=max Fla.r)—f(r)\. is odd or even) at a distance less ‘than e for jany, 

, e>0. Clearly this osmmaonaahion may be modified 

so is F(a,r). This contradicts the validity of theorem | to obtain a point in 7” on the &th axis with opposite 
B. i sign. Hence 7’ is axial. 
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Lemma 3: Let P’ be a closed arial set; then P’=E,,. 

Proor: The proof is based on the following as- 
sertion: Every closed n-sphere whose boundary passes 
through the origin in /,,; contains a segment of a 
coordinate axis which is attached to the origin. 
That this assertion is true is seen as follows. Con- 
sider the tangent plane of the sphere at the origin. 
If the positive #th axis is on one side of the plane, 
then the negative Ath axis is on the other side; i.e., 
at, >0 implies a,(—1,) 0. Now the tangent 
plane cannot contain every axis; hence there is an 
axis intersecting the tangent plane at the origin. 
Clearly this axis must intersect the sphere also, which 
establishes the assertion. 

If P’ is closed, the complement of P” is open. 
Let ¢ be a point in the complement of 7’, and let S 
be the smallest closed sphere with center at ¢ which 
intersects 2’, and let 6 be a point of intersection. 
Now S contains a segment of some coordinate axis 
of a coordinate system with origin at 6. Since 7’ 
is axial, that segment contains a point of P?” which 
contradicts the definition of S and hence the ex- 
istence of ¢. 

With these three lemmas the second question may 
be answered by the following theorem. 

THEOREM 2: Theorems A and B are both valid for F 
ifand only if P=E,, and Tis a Tchebycheff set. 
~ Proor: It is a classical result that theorems A 
and B are valid if P=F, and TJ 1s a Tchebycheff 


set. 
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From lemma 1 it is seen that the validity of 
theorem B implies that F has property Z. If F 
has property Z and theorem A is valid, then theorem 
| implies that ? is closed. Now there is a one-to-one 
mapping between P? and P’, the set defined for 
lemma 2, which is continuous both ways, since / 
has property Z. Hence, if ? is closed, then so is P’. 
From lemma 2 it follows that P?’ is an axial set, 
and lemma 3 implies that ?’=F,. Hence P= F, and 
T is a Tchebycheff set. 

4. The answer to the third question is given by 
theorem 3. The proof of this theorem ts rather 
complicated and since it is aspecial case of a theorem 
to be published elsewhere, the proof will not be 
given here. The result is given here for the sake 
of completeness. 

Tueorem 3: Theorem B is valid for F if and only 
if (1) F has property Z; 

(27) given a*éeP, k< n, 1 ,\0=2,.< 7; =, Gonna 
Teyy= 1} and € with Oe 4 min (4)4,—2;), = 9, 1, 
, hk, then (a) there are a,,a,€P such that, for 
(0,1), F(a*, 4) —e< F(a,,24)< F(a* x)< F(ay,2) 
F(a* x) +, (b) there are a3,a,€ P such that | F(a3,2)— 
F(a*.2)\<\e, |F(ayr)—F(a*,x)\ Ce, for x €|0,1) and 
F(a3,2)—F(a*,x), F(ay,2)—F(a*,x) change sign from 
xr,—e to 7;+€, } i. ie « .. k and have no zeros 
outside |x,—e, 2)+€|. Further F(a3,0)>F(a*,0)> 
F’(a,,0). 


te 


Wasuincton, D.C. (Paper 64B2-25) 
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Note on the Solution of Riccati’s Differential Equation 
H. Herbert Howe 


(February 4, 


1960) 


Three recurrence formulas are developed, giving the solution of a particular case of 


Riceati’s equation in power series valid in the neighborhoods of 0, ©, 
The first two were programed for computation on the SEAC, 


point , respect i\ ely. 


and an arbitrary 


1. Introduction 


A particular case of Riceati’s differential equation 
is important in the theory of radio propagation. 
In the notation often used in that subject, omitting 
subseripts irrelevant in the present paper, the 
equation Is 


a oy, 41=0 
dr 


(1) 


where 6 will be considered the independent variable, 
and both variables are complex. The equation 
may be written 
dr l = 
dé 26°r—1 


Since this is a first-order differential equation, 
there should be one constant of integration. — If 
the value of 7 is arbitrarily chosen for some value 
of 6, the differential equation ought to define + for 
all other values of 6. We may, for example, arbi- 
trarily select 7(0) or 7(@), and from it derive + for 
all other points. 

If there exists in the complex 6-plane any point 
P such that as we approach it the value of 67 
approaches 3, then the derivative dr/dé approaches 
infinity, and 7? is a singular point of 7. Since 67 
involves both the independent and the dependent 
variables, the position of ? depends upon the par- 
ticular function 7+ which we are studying: e.g., 
upon the initially selected value 7(0). It not 
immediately obvious that any such /? exists, nor 
that there may not be more than one. 

A single-valued function of a complex variable 4, 
with isolated singular points, can be represented 
[1] by power series as follows: In the 6-plane, draw 
concentric circles, one through each singular point, 
and all centered at the origin. Within the smallest 
circle, the funetion is represented by an ascending 
series in 6 (i.e., there is a constant term, and a series 
of all positive integral powers of 6). This series is 


Is 


absolutely convergent and represents the given 
function for all values of 6 within the circle. Outside 


the smallest circle, it diverges. On the circle itself, 
it may or may not converge. This circle is called 
the “circle of convergence” for that 
Special case, its radius may be zero. 


series; as a 
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Outside the largest cirele, the function is repre- 
sented by a descending series, namely a series con- 
sisting of a constant and negative integral powers 
of 6. It is absolutely convergent outside the circle, 
divergent inside it, and may or may not converge 
on it. 

In the region between any pair of circles, the 
function is represented by a series containing all 
integral powers of 6, both positive and negative, 
and a different series is required for each inter- 
circle region. Such a series is absolutely convergent 
in the region between the two circles to which it 
appertains, may or may not converge on those 
circles, and diverges outside the specified region. 

All of these series together represent the given 
function at all points of the complex plane, except 
perhaps those on the various circles. The function 
exists for all points except the singular points 
themselves; for other points on the various circles, 
one of the series may or may not represent the 
function. 

For the function in which we are interested, we at 
first use the hypothesis that there is one singular 
point, and that therefore there are two series: The 
ascending and descending ones, which have a com- 
mon circle of convergence. A few terms of each series 
are given by Bremmer [2], and a greater number of 
terms are given by Johler, Walters, and Lilley [3]. 

Although one of the two series converges for any 
value of 6 that does not lie on the common circle of 
convergence, the rate of convergence is slower the 
nearer the point is to that cirele. For points very 
near the circle, the converging series will not, in 
general, provide a practical method of computation. 
Even within the limits of practical computation, it 
will be impossible to select any fixed number of terms 
at which to stop the computation, if given accuracy 
of the computed result is desired. Rather, the num- 
ber of terms must be increased as 6 gets nearer the 
circle. 

One could, of course, always compute a fixed num- 
ber of terms; but there would then be no assurance as 
to how accurate the result would be, nor as to which 
series should be used, nor even that either one con- 
verged. There is thus need for a formulation of the 
series which does not rely on expressions with a fixed 
number of terms. 








2. Ascending Series 


In 1955, the present author programed radio appli- 
cations of (1) on the National Bureau of Standards 
Electronic Automatic Computer (SEAC). Some of 
the results have been published [4], but not the for- 
mulas. Instead of using a fixed number of terms of 
the series, recurrence formulas were developed for 
expressing each term as a quadratic function of the 


preceding ones. To derive the formulas for the 
ascending series, We express 7 as a series in 6: 
r=») 5,8". (3) 


n=0 
Rewrite Riceati’s equation as follows: 


| 


ro (4) 


26°r+-1=0. 


Insert into (4) the value of 7 given in (3), and in one 
term change the summation index from n to m: 


—2 S) 5b,,6"*?7+-1=—0. 


| m=0 (>) 


Now clear of fractions, remembering that two abso- 
lutely convergent series may be multiplied term by 
term: 


(). ()) 


1-2 >) nb, bd" "14 nb, 8" 


m=0 n=0 0 


In the second term, introduce a new summation 
index k, and eliminate the index m, & being defined by 
m=k—n. To find the summation limits for n and & 
make a plot of m and n; in (6) we are to evaluate the 
term for each point of the following grid: 


© ;) 
] 4, 
5. 
9) 
* 

0). , . ; 

0 ] 2 ) 4 7) 6 

Wi x 


Since k=m-+-n, lines of equal & are the diagonal 
lines. Since we need to sum for all points of the 
grid, extending upward and to the right indefinitely, 
it is seen that we get the same points if we let & run 
from 0 to , and for each value of & let » run from 
0 to k. 


> k fon) 
1-259 >) nb bend’ +5 nb, 5"'=0. (7) 
k=0 2n=0 n=0 
Now change the summation indices again. In 


term replace n by /. 
terms for ) 
common summation: 


Then write out explicitly the 


(S) 


1+ b+r[ J) 2> 5 nb,b,_, Jo 

Since this equation is to be satisfied for every 

value of 6 within the circle of convergence, the 

coefficient of each power of 6 must vanish. This 

imposes no condition on 6,, which is equal to 7(Q). 

which we have previously proposed as the constant 
of integration. From the constant term of (8), 


b, I (9) 
and for the higher powers of 6, we that the 


bracket of (8) must be zero for every value of 7>2, 
L.e., 


See 


(10) 


j-2 
h (2/7)> 5 nb,b 


In particular, b,—0. 
Now break the right side of (10) into two equal 
sums, each having a coefficient of 1/7; and in one 


of them replace nN by a Hew summation index m 

defined by n=j-2-—m: 

, 1&2 =a 

) > > nb b aT: 2, (7-2 m)b, 5 »b,. (1) 
J n=0 Jom 2 


Now replace summation-index MW by NM. Since the 
summations run over the same values of n, and since 
terms with the same » have also the same b-product, 
the series can be combined to give: 


‘ a 9 
h J ; — ; bb 
J n=0 


» for 722. (12) 


Equations (9) and (12) show the coefficients to be | 


used in (3), after the summation-index of the latter 
has been changed to }. 


3. Recurrence Formulas for the Terms 


For the SEAC computation, it is better to use 
recurrence formulas involy ing the terms rather than 


the coefficients. Let 4; represent a term of our 
series: 1.e., 
B h 6’. (13) 
Substituting into (3), (9), and (12), we get 
‘ | 
r= >, BB, B, 6: B,—0: 
0 (14) 
7—2 .S pp ,; ' 
B anes i > BoB for J2>, 
} sues | 
where By=—b)—7(0) is the constant of integration, 


arbitrarily. In practice, it is a 


and may be chosen 
chosen to fit the conditions of 


boundary condition, 


the second term, replace k: by } 2, and in the last | the problem at hand. 
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ly the 4. Descending Series 
‘ing my 
7 Using a similar but somewhat more complicated 
method, we get for the descending series 
(X) | | ‘5 
> es ‘ 
; 244A Ai ( a4,)\3): 
every 
, the l [—! e 
~—: h eal » 15 
his Ay=| 2A,) Boy Ants ,: 
T(Q), =e 
iStant  . — onl | for jz 
; J 0 a 
(9) where Ap=7(@) is the constant of integration. 
In order to determine which value of sy belongs 
t the’ t the same function 7(6) as does a given value of 
J2 B,, so that the results of the two series may be com- 
hined into a single function, other information about 
(19) the problem at hand is used [3, pp. 1 to 4}. 


5. The SEAC Computations 


5 ‘ , ‘ - . 
equal The SEAC code used (14) and (15) as follows: 


1 one (ne of them was selected more or less arbitrarily. 
eX m The terms -1, or B; were computed and saved, since 
each involves those which come before. When 
seaside, computations were continued until a term 
(yj) was found which was negligibly small. (Real and 
imaginary parts of 1, or , each smaller than, say, 
' 10°.) This determination began with j=3, since 
» the B, is identically zero. When such a term was 
since reached, the computations were discontinued, and 
luet, the A’s or B's were summed. 
Ifany .1, or #, (starting with v1, or By) was very 
large, that computation was abandoned, and the 
19) other formula was tried. The same thing was done 
7 if no negligibly small term was found within the 
space available for storing terms in the memory 
obe | (about 65 terms could be stored). If, by the same 
utter eriteria, the second formula did not converge, the 
computation for that particular value of 7(0) was 
abandoned. Some of the gaps in the published 
; results [4! come from this cause. (Some occurred 
because of lack of time to complete the lengthy 
use COMputations that had been planned.) This might 
than occur either because the point was very near the circle 
our of convergence, so that still more terms were needed; 
or because for certain functions there are more than 
(13) one singular point, giving a region in which both 
series actually diverge. 
6. Convergence of the Series 
The criterion used for terminating the computa- 
(14) tions (real and imaginary parts of some term both 
less than 10°°) is not rigorous, and under some 
circumstances it would be incorrect. This may be 
most easily studied by substituting into (12) so 
ion, | that each 6, is a polynomial in bo= Bo, to give the 
Is : » equations of [3]. If By happens to be such that any 
s 0 
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b, is zero, then B, will be zero regardless of 6, even 
though subsequent values of B might be of sub- 
stantial size. For example, two out of every three 
b’s are zero if By=0; or 6, is zero if 54+8B3—0, ete. 

For the problem under consideration, each value 
of By for which 7(6) was to be derived had an ampli- 
tude of 7/3, and the moduli were all transcendental 
numbers exceeding 1.8 [3, p. 20]. Being transcen- 
dental, they could not make any 6, exactly zero; 
and it is extremely unlikely that, among the small 
number of values that were used, there would be 
one that would make any 6, negligibly small. The 
expansions of 6, through j=11 as polynomials in 
by |3| show that ‘the coefficients of powers of bp tend 
to increase as j increases; and since the modulus of 
by exceeds unity, it seems most probable that the 


values of 6, continue to increase as j increases. 
Hence, it seems warranted to conclude that when 
any 6, becomes negligibly small it is probably 
because 6) is less than unity, and {|6/| is small; so 


that succeeding values of |B; probably will be still 
smaller, on account of the factor 6 in (14). 

A similar argument seems to apply to the descend- 
ing series. More reliable criteria would be to require 
that each of several consecutive terms of the series 
be negligibly small. 


7. Bridging the Gap 


Although neither series can be used for practical 
computation for points near the circle of convergence 
(or between the circles of convergence, if there is 
more than one singular point), the function itself 


is regular in any region not including a singular 
point; and for a point far from a singular point, 


the derivatives would not be excessive. Hence, if 
we are far from a singular point, it may be entirely 
proper to interpolate. across the gap where neither 
series could be used. 

Another method getting values in this gap is 
by means of analytic continuation. We can use the 
ascending series to give a value of 7 for some value 
of 6 which we may call It is then possible [1, 
p. 196] to represent 7 as an ascending series in 6—q, 
which will be valid within the circle centered at gq 


and passing through the nearest singular point. 
Such a series, obtained by methods analogous to 


those already deseribed, is 





“ , : 5 —¢ 7 
r= > (o=7(Q); G : 4 : 
j=0 200g —1 
) 9 al 2 
At 2 0q°) . (6 qg)? >> 6,0; n—-2 (16) 
n=0 
sf) eee ' ji—1 
2q(6—q) / F »p & CC j-e-17 >> C,C;-» 
Jj n=0 n=1 Z 
No computations have actually been made with 


this formula. 
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+ band On a Generalization of the Index Notation for 
Absolute Tensors of Arbitrary Order 


re. 574 
2-96) | Edmund H. Brown 


(January 7, 1960) 


A generalized vector index notation is introduced which facilitates study of properties 
of tensors regardless of order and variance of components. This generality is obtained by 
replacing the set of m indices, each of which vary in a single domain of n integers, by one 
vector index which varies over the integral lattice points of an m-dimensional domain. The 
notation also suggests an elegant treatment of order transformations of tensors. 


1. Introduction 


The usual index notation for the components of a tensor in a given coordinate system has 
the advantage that repetitious display of similar components of tensor quantities or equations 
is avoided without sacrificing directness and simplicity of algebraic manipulations. At the 
same time, retention of a single base symbol (to which the indices are attached) for each tensor 
is conceptually appealing since there is a one-to-one correspondence between base symbols and 
covariant objects. This correspondence is especially desirable in the expression of physical 
theories since physical observables are required to have the covariance property. A disadvan- 
tage of the notation is that, in the study of general properties of tensors, equations such as the 
transformation relations, the expressions for covariant derivatives, or the integrability condi- 
tions do not retain the same form, but require an increasing number of terms or factors with 
increasing tensor order. In this note, an amusing generalization of the usual notation is pre- 
sented which permits a study of such general tensor properties in a form which is invariant 
with respect to order. In addition, this generalization suggests some new concepts, such as 


that of the order transformation. 
2. Index Notation 


Let 7 be an mth order tensor defined on an n-dimensional coordinate-space G and z' 
be an arbitrary coordinate system in G. Then, m indices j, may be defined in the domain 


J(,) consisting of the set of positive integers such that 0<j,<n fora=1,2,. . .m. By assigning 
one component of 7 in the coordinate system sz! to each set of the indices j, the tensor 7’ may 


be completely specified and written in the index notation as [412 >>» 4m 
For every coordinate transformation in G@ given by a set of n analytic functions z°=<s"(zr‘) 


there will be a transformation of the components of the tensor 7 at a point 2‘ given by 


g8182 ++ %m— piila---ImyMip%? | |, zr”, (1.m) 
where 
s, Or*e 
I * 
¢ 6Or" 
Thus. for a vector, m=—1, the transformation law will be 
t*(2r’) t?(x*)z3; (1.1) 


and, in general, the transformation law of an mth order tensor will involve exactly m of the 
. . 1 
transformation matrix factors s,. 
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If a tensor field 7(G) is defined on all points of G, the true rate of change of Tis given by 
the covariant derivative 
Seite h 
{772 = Sr +t" , rat eh rs? (2.m) 
where the convenient vertical bar notation of Green and Zerna' is used. For a vector. the 


expression for the covariant derivative becomes merely 


Of! 


at OT nn’ (2.1 ) 


and, in general, the covariant derivative of an mth order tensor will involve exactly m “correec- 
tive’? terms in the Christoffel Connection I‘). 
As a final example, the integrability conditions for an mth order tensor may be written 


i i,i i hi i » i,h ry } h [im ‘ 
t! | ee a m t" mi" att ’ mA? gles. HO ie. (3.m) 


where /'),, is the Riemann-Christoffel curvature tensor. For a veetor, the integrability 


condition reduces to 
t*|, —t'| gO R yn: (3.1) 


and, in general, there will be exactly m terms on the right side of the integrability conditions 


(3.m) for every mth order tensor. 
3. Index Vector 


In the previous section, exactly m indices /, were adjoined to the base symbol of an mth 
order tensor, with each of the ), varying in a single domain J,,).. We now consider the converse 
possibility of a Cartesian lattice index space J,,,, which is the mth Cartesian product of m of 
the spaces J,,, and define a Cartesian indes vector jin this space. For the moment, we consider 
only index vectors whose components are positive integers, that is, we restrict J,,,, to the first 
hyperquadrant. For an mth order tensor on an n-dimensional space, a marimum index vector 
jcan be defined by the relations j,=", a=1, ... , m, for its components j,. Then, since 
there will be exactly n” index vectors j which satisfy the conditions 0~— /j,</,, one component 
of the tensor may be assigned to each of index vectors j and the tensor itsell represented by the 


symbol ?’. 
One of the simplest methods of forming higher-order tensors is by adjunction of vectors. 
Since coordinate di the rentials are vectors, a first application of this index vector notation results 


in the symbol da’=da'idr? . 2. dz’. Although the transformation matrices +, do not repre- 
sent tensors, the svmbol sj=sja 2. . 2," may be constructed by analogy with da’. Using 


this generalized matrix, the transformation law for a tensor of arbitrary order may be written 


Arf, (1) 
where (corresponding to the usual summation convention) the diagonal repetition of j indicates 
summation over all the index vectors in J, 

By analogy with the derivation of Brillouin,’ generalized expressions for the Christoffel 
Connection may be obtained by considering the displacement of a tensor of arbitrary order. 
If 2 are normal coordinates, that is, coordinates such that the covariant differential D#* is 
equal to the ordinary differential df in the infinitesimal neighborhood of the point «?, then 


t*(7°+-dr’)—t*(r")4 Dr 


A. E. Green and W. Zerna, Theoretical elasticity (Oxford University Press, Oxford, England, 194 
2 L. Brillouin, Les tenseurs en mécanique et en élasticité (Dover Publications, New York, N.Y., 14 
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Transforming to arbitrary coordinates 2' by means of the transformation matrices z,(2°+ dr") 


and s(x”) gives 
tirt)—t8(2")aa(2"),  O(a*® +-dr*) = [t°(2") + Dt") h(a? +d’) 


f(r") +Dt* Ta(r”) + Taylr” t'(2*)+-Dt'4 triprsuida", 
where ai, =0.4/0r? and where we neglect the term containing the second-order differential, 
Diry,dz’. Defining the Christoffel Connection by Ty = — ryrert, the covariant differential 
can be written 


Dt’ dt’ T PT da", 
and the generalization of the expression (2.m) for the covariant derivative of a tensor of 
arbitrary order becomes 
: or — 
i Jy ») 
t |, oer! rT (2) 


If the normal coordinates 2* are assumed to be obtained by transformation from a third 
set of arbitrary coordinates 2’, we can write 


‘i i >t ay? © ,ap? pi,;p,p* A apt pay? p? 
ri (reer Pipl ots tel PLk— ST, Teter rk 
or, 
‘i 1° pi pPpt — pli pep? 
, rhr?. ea 
hy Ns Hae eel the 


which is the generalization of the well-known Christoffel transformation law. Clearly, the 
original definition of 2} requires that 


Ar, / Q--° ‘mph 4 / Ivers m|" ine 2 6an f*a"gees [tm 


“jk “jks 


ia 


by the law of differentiation of products and the relation 6)=.242%. 
For an arbitrary tensor field t'(r*), the expression for the covariant differential permits 


a calculation of the change in the components at a fixed point after integration around a given 
curve Cof G. If é6f is the change at the point, then 


éti p dti p Dir p phy dx ( ti\,.da*— poride 


( c 


( ea (/ Om, OT 
2J. Ou ou" 2J. ou’ Ou" 
l¢ a : : wis ‘Ort, THe a) ' ‘or, “om . 


where S isjthe surface enclosed by C and ds’* the differential surface element. Following the 
usual definition ofthe’curvature tensor, 


ri or! om ORs owe a 
Py) k ( —+ 4 in) —ere “i ty) 


so that the change 6f may be written 


— 
- 
- 


a || { til p— tha, +t Rig jds®, 








Since integrability requires 6f'=0 for arbitrary area S, we arrive at the generalization of 


the integrability conditions (3.m) in the form 


7, —tp=—OR (3) 
Thus, to be integrable, a tensor field cannot be completely arbitrary, but must satisfy the con- 
ditions (3). In the case of a tensor field obtained by parallel displacement to all points of a 
tensor given initially at one point, we must have )f'=0; thus, the field can be integrable only 


if the space is Euclidean, Ri, zp =O0. 


4. Mixed Index Vectors 


The index vector notation facilitates study of general tensor properties regardless of order. 
In such a study, however, specification of the covariancy or contravariancy of components with 
respect to particular indices may be even less important, especially in a metric space where 
components may be immediately transformed from one type of variance to the other. Thus, 
a general notation should allow interpretation of such differences in variance without actually 
displaying the differences in its symbols. 

Let us consider what significance can be attached to an index space -J,,,., which includes all 
hyperquadrants, that is, which allows both positive and negative integers as components of the 
index vector j. Each hyperquadrant in the enlarged ./J,,,, will be bounded by m hyperplanes 
formed on m semi-axes and, thus, there will be 2” such hyperquadrants in J,,,. Defining a 
maximum index vector j by j,=-En with some particular signature o , +, 
for the signs of the components then, with the restrictions sign j,=sign jaand 0< |je <\j4|, there 
will be exactly n” index vectors j for each maximum index vector], that is, exactly the same 
number as components of an mth order tensor defined on an n-dimensional space. In addition, 
the number of maximum index vectors will be, clearly, equal to the total number of combina- 
tions of covariant and contravariant indices for an mth order tensor. 

The above equality, thus, permits the following correspondence: To each contravariant 
index j, we assign a positive sign to the respective component of j; to each covariant index j» 
we assign a negative sign to the respective component of j; and, as the index vector j varies in 
Jim, subject to the conditions sign j,—sign Ja and 0< jai < Ja , We assign to each j one of the 
respective components of the mixed tensor corresponding to j. All the 27" index vectors in Jin) 
thus specify all the possible covariant, contravariant, and mixed components of a given mth 
order tensor defined on an n-dimensional coordinate space. 

In particular, if # are components of arbitrary variance (that is, some of the components 
of the index vector may be negative) of a given tensor 7, then t '=4, and the form of the trans- 
formation law (1) remains valid for this generalization if we also define 2f=.t. Thus, (1) 


may be written as either f=fzj or t, try regardless of the variance of the components. 


el 


k st Stn , . ’ 
Iyj— 9, 8. . . O™, where the 


> } ij tJ i } lm m<¢ si 
Defining dij 9i,j Ii,j, — Jij. I” d oe , and 0; q 


indices 7, and j, may be either positive or negative, the corresponding ‘“‘covariant’’ Christoffel 
° ; . ch - , , 
Connection will be given by Ty jx =gin jx. Using the identity 


OO”; 
“ik » ‘ | . | | " 
or (Tye, Dy ji, Gir, ‘2a Di ken eo“ (Ty ix, 4 ys in) Vii sd So qj j si 
the expression (2) for the covariant derivative of a tensor of arbitrary order may be written 
in either of the forms 


r or" Pwr , Of; sk 
: arte Ce ol tii a Se 


regardless of variance. The similar generalization of the integrability conditions (3) for ten- 
sors of arbitrary mixed variance is obvious. 
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5. Order Transformations 


The identity gj —t-' which follows from the results of the last section suggests the definition 
(t)-'—t-' or, generalizing the latter, the definition f'=(#)*=t9, where @ is a rotation operator 
in J,») Subject to the restriction O— |i,;Sn. Thus, a variance transformation of the compo- 
nents of a tensor 7’ on G can be considered equivalent to a restricted orthogonal transformation 
of j in the index space -J,,,). 

By removing the restriction, a more general orthogonal transformation @ can be defined 
which can be called a decreasing order transformation, that is, @ will be a transformation which 
rotates j into a subspace J,,) of J:,) and, correspondingly, transforms the tensor 7',,) into 
7... If a@isa marimum order transformation, that is, one which transforms j into a zero vector 
at the origin of J,,,), then clearly there will be the same total number of such possible maxi- 
mum order transformations @,,, as there are invariants of the tensor 7. 

In particular, the sealar f=fit-'=f4—te'--t;,...;, is the only invariant possessed 
by tensors of a// orders (and thus might be called the ‘fundamental quadratic invariant’’). 
Just as contraction of a tensor is indicated in the ordinary notation by summation over a pair 
of indices (one covariant, one contravariant), contraction in this new notation results in a 
tensor fie =X), ifm, subject to the restriction ¢,+7,—0, with k, h fixed, where v=n—2. 
The reduction in dimension of the index space results from, first, restricting the index vector 
to the hypersurface defined by 7,+-7,=0 and, second, summing the components of f' (and as- 
signing each sum to the corresponding sum of the index vectors) over the values of the index 7,; 
thus, the dimension of the index space is reduced by lwo. 

Further consideration of general order transformations should be of considerable value in 
the theory of groups as used in quantum mechanics, in the study of irreducible tensors, and in 
determining the invariants of tensors of arbitrary order. If such general order transforma- 
tions were unrestricted they would clearly lead to tensors of more general type, that is, tensors 
represented by a set of components defined, not as functions of a point in a single space, but as 
functions of VW points in WZ spaces having possibly different dimensions. 


Bou.per, Coo. (Paper 64B2-27) 
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Upper and Lower Bounds for the Center of Flexure 


Lawrence E. Payne 
(February 1, 1960) 
There is disagreement in the literature as to how the center of flexure of an isotropic 


elastic beam should be defined. In either of the two most widely accepted definitions, upper 
and lower bounds for the coordinates of the center of flexure are obtained. 


1. Introduction 


We consider a cantilever beam of uniform cross section and choose the z-axis to coincide 


with the line of cross-section centroids. The plane end 2=0 is assumed fixed and the other 
end 2=/ is loaded by some distribution of forces which is statically equivalent to a single force 


with components (W,, W,, 0) acting at a load point (7,, y,,/). Inthe treatment of the flexure 
problem |9|! (a similar treatment appears in [3]) we seek a load point (z,, y,, /) which has the 
property that any load (W,, W,, 0) applied through this point produces no local twist at the 
centroid of the section, i.e., the mean value of the local twist over the cross section vanishes 
(see [9] p. 200). This point (zy, y,, /) is referred to as the center of flexure [9]. The determina- 
tion of (vy, yy, /) permits us to break up the general flexure problem into two separate problems 
one of pure torsion and one of pure flexure. This point is sometimes referred to as the center 
of shear (see, for instance, [3]), and it is well known that there is difference of opinion in the 
literature as to the way in which the point should be defined (see, for instance, [1, 2,4, 10, 11). 
A discussion of these differences has been given in a recent paper by Pearson [7], who demon- 
strates that in general the two expressions for the center of shear vield results which are very 
near to one another. 

It is not our point here to enter into the discussion of the relative merits of the two defini- 
tions. We merely point out that the expression of Trefftz [10] is simpler and independent of 
Poisson's ratio. The center of flexure of [9], on the other hand, permits an easy decomposition 
of the general flexure problem. 

In order to avoid confusion we shall hereafter refer to the point defined in [9| as the center 
as the center of shear. 





of flexure and that given by Trefftz [10 
In this paper then we derive upper and lower bounds for the coordinates of the center of 

flexure (ry, yy, /) and the center of shear (x,, y,,/). In the first case we use the formulation of the 

flexure problem derived in |9| and thus concern ourselves with the problem of pure flexure. 

It is of course well known (see, for instance, {[7]) that if we have the solution (o the pure 
torsion problem for the beam in question then we can determine both the center of shear and 
the center of flexure without knowledge of the solution to the flexure problem. As would be 
expected it is possible in each case to obtain upper and lower bounds for the coordinates of the 
point by approximating the solution to the torsion problem alone. 


2. Bounds for the Center of Flexure 


We denote by 4% the cross section of the beam and let C be its boundary. (The region 
may be simply or multiply connected.) The origin of a rectangular coordinate system is 
taken at the centroid of Z. For simplicity the 2 and y axes are chosen to be principal axes of 


inertia, Le., 
. . 
ryd A=0. 
: 9 
I (2.1) 
G 
1 Figures in brackets indicate the literature references at the end of this paper. 
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According to [9] the pure flexure problem is solved once the solutions to the following boundary 


value problems have been determined: 
Oy) Ory; _ Oy, , Oe» 


s ,=0, ===() in ¢ 
Or OY? Or" T >) — 
oe [((1+o)s2*?—oy*|n, 
[((l+o)y?—or?|n, on OC 
‘ 


where ¢ is Poisson’s ratio and n, and n, are the # and y components of the unit normal on ©. 


(The normal is assumed directed outward from %.) The coordinates of the center of flexure 


(r,, yy) are defined by: 


. Po (oom o = 
vy =xTTeT. IJ Ov —_— 7. —(l+o)ry or \dA, 
dv, dv, 
= | ae | 2 sl, 
a “(1 } on dso 1 Ov y or Hl +o)e'y—oy"d A 
where 
| | yu 
‘GY 
| |. “AA. 
vy 
Thus 


Using the divergence theorem we have 


2(1+o)/,r,= Pees (rny—yne)ds J 


( 


where 
. [((1+-o) ry?—or? |\d A. 
7 


Since J is a geometric quantity it can be computed explicitly. 


2(1+o)L,7, \J | 2 (rg.)— 2 (ye2) — (+a) ry? os | LA. 
Y 


(2.9) 


(2.10) 


We note now that the coefficient of g in the boundary integral in (2.9) is proportional 
to the boundary value of the normal derivative of the warping function ¢ in the pure torsion 


problem. In fact 
: yi n ( 
U Th, 0 e 
yn Y r y 


Thus if we let 


x=2(1+06)/,r,4+-J 


-> ¢2 oF ds —D)(¢,¢2) 


we find 


where 


ite a) Oy Ov, +50 eV dt, 


Or Or Ty OY 


(2.11) 


(2.12) 


(2.14) 


The determination of bounds for zs, is thus equivalent to the determination of bounds for x 
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Note that if the warping function ¢ is known the value of z, can be determined explicitly 
since 
aoe | gl(1+-0)y?—or*|ndls. (2.15) 
c 


If ¢ is not known we introduce two arbitrary harmonic functions vu and uw, and define 


a,=—D(¢,uz) —Dus(yn, rn, )ds, (2.16) 
C 

a,.=—D(u,@) — ful(1+-oy? ox?|n,ds, 2.17) 
Cc 

d3—=—D(u,uy)- (2.18) 


», and a, explicitly. We now form the expression 


With wand ww prescribed we can determine ay, a 


X—a,—a,+a, D(g—U,go— Ue) 


A Ou» ‘ 
—4 (~ 0{ (0 +-o)y?—ox"|n, on pus (2.19) 
3 


An application of Schwarz’s inequality vields 
— 4 : ' Ou, \ ° 9 « 
(X—, —d)+ a3)? << De ujtdsp \(1+-0)y- or’ |ny— > ds. (2.20) 
n 
GC Cc 


It remains now to approximate the first boundary integral on the right of (2.20) in terms of 
known quantities. 

By adding an appropriate constant (if necessary) it is possible to choose g-« in such a 
a way that 


Pie —u)ds=0. (2.21) 
Cc 


With ¢—u so chosen we shall estimate the first integral on the right of (2.20) in terms of the 
boundary integral of the square of its normal derivative. To this end we note that 


P(e) f(Z)u 


a > min “—, 
fp (g—u)*ds Pris =0 Prids 
“ 


. 


A 
where L denotes the length of the perimeter of C, 0/0s denotes the tangential derivative on 


(’, and v is any function continuous on ( and satisfying the condition p rds=0. The quantity 
on the right of (2.22) is just the eigenvalue of the vibrating string equation, the differential 


equation which arises as the Euler equation for the minimum of the quotient (Pia ds) ds ) 


C 


Peds. It follows then that 


o 
eee 2 | 
Pie u)'d8 <7 p | °. (y- | ds. (2.23) 
Cc Cc 


We seek now an inequality which relates the integral on the right-hand side of (2.23) to the 
boundary integral of the square of the normal derivative of g—u. We introduce therefore, 
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an arbitrary continuous vector function with components f/, and consider the generalized 
Rellich identity (see [6, 8]). Letting w=¢g—u, we have 


*(. Ow ~ . Ow ow S "Of; ow oY “Of, 
{fp Sz, Awd. A | th <= ds—; Sere grad w “_ Jar ‘ 3) | . grad w\dA 
. “ (2.24) 


where A denotes the Laplace operator (A= (0*/0r*)+(07/O0y*?)) and summation is to be 
carried out over the repeated indices k,/=1,2. The f, are now assumed to be so chosen that 

p Any -Oon C. (2.25) 
A decomposition of the derivatives of uw in the boundary integrals into normal and tangential 


components vields (see [S]) 


P % ow \* 4 ow? a * Ov, OW OW Col, = bow oe 
| p( = ) ds Po - ) ds 2q f, oo a ds+ {J se grad wd. ‘J | 3 dA. 


Cc C c 
(2.26) 
Since the functions f, are prescribed, we can find a constant y such that (see [6, p. 595]) 
"COS, ’ "Of, r+ oe 
| | : grad wd. *IJor , das yD ww (2.97) 
We make use also of the ordinary Green’s identity 
~ ow 
Dwyw) | w ds, (2. 98) 
On 


which vields by Schwarz’s inequality 


2 1 


D(w,w) <( | ur ds y( p( e. ) ds y’ (2.99) 


- 


( S 


In view of (2.23), eq (2.29) vields 
L (=) > 
D(ww) <5 LY > (2.30) 
<7 i] 


We now insert (2.27) and (2.30) back into (2.26), make use of the Schwarz inequality and 


obtain 
~ (Ow? ~ (ow? | ( S ) ow \? ~ (ow? , 
Pr( ds )as<f v( on) +2 f p (3) L Pol Os Jan] 
‘ks Or — 
tO epi Maal ; ol if ( on) as} ar 


Where Pin is the minimum value of pon CC. This inequality vields 


OK \ ; 1 Or, \ , 
( ‘ e) 1 , { l S ) ow ds 5 i ’ 2 1 | | S ) w\- ; 
lf o( ep (on) +L ara (Fy LP an) 
} , 


p 
2 : 1 
' +4 “Le? (=) )as )+ ‘2 Dp yds *, (2.52) 
4a” Pi ns (on on 
3 
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From (2.32) it is thus always possible to obtain the inequality 


nd ow ‘ “ye 
p( Os ) ds <a f ( on ) ds (2.33) 
( ri 
where 
@Pmin + [B°+ Pmax]? 3)? (2.34) 
and 
,Or) 
Os 
P + (2.35) 
p* max 2 Pmin 


If the boundary ( is star-shaped with respect to some point (z,,y.) in %, the problem is 
simplified considerably. The fact that the boundary is star-shaped with respect to (#,,y,) is 
equivalent to the condition that 


q=(4—29) Ne + (Y—Yo)n, >0 (2.36) 
at every point on C. Thus, in this case, (2.25) is satisfied if we let 
fi=7—--d, fo=Y—Yo. (2.37) 


The left-hand side of (2.27) then vanishes identically, and hence we take y=0 in (2.31), (2.32), 


and (2.35). We find in this case instead of (2.33) (see [5]) that 


b (So) aes aaa [ 1+ 3 | | $5 (2°). ds. ( 


In the general case then, an inequality of the form 


~ (Ap—uV, - £ e) ya 9 26 
( ds ) ds < A E (y » | ds (2.39) 
ps 


is obtained, where 1 is a completely determined positive function. Inserting (2.39) and 
(2.23) into (2.20), we find 


 ' & Ou 2 ~ f Ou, _ . ‘ 
(X—a,—d.+d3)*< te ¢ A | Syn. 7H ] ds t { ame to)ye cx'n, } ds. 2.40) 
( rf 


(It is apparent that the quantity 4 may be placed in the second integral rather than in the 
first if it is desirable.) We now use the Rayleigh-Ritz technique to make the two integrals 
on the right of (2.40) small. 

Note that it is possible to derive close upper and lower bounds for 2, without approximating 
the flexure function at all. In fact, if we use instead of ua conjugate function v defined by? 


to 
3 


( 


tf 


ou oO Ou ov (2.41) 


or ~ Oy’ OY Or 
we obtain (setting u.=0) 
- a fee . ° ' 
X— @\° < : iD (7 4?) ds 4) AA 11 o)y ox*|n ly “ds. (2.42) 
4n* - Os h 
c C 
In this case the .1 term has been put in the second integral since this expression is a purely 
geometrical quantity (independent of 7) and can be computed explicitly. 
lhe function p is not to be confused with that used in (2.22). 
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In an entirely analogous way we find that 


2(1+o6)L,y,4 (1+ o)2°y—oy' |\dA=D(¢,¢)). (2.43) 
Yo: 
By introducing the two arbitrary harmonic functions uw and wv, and proceeding as before, 
we obtain upper and lower bounds for y,. Note that if we choose w,,7,=0 and merely approxi- 


mate the warping function (or the conjugate function) on the right, then the same minimizing 
function uw (or 7) minimizes the error term in each case. 


3. Pounds for the Center of Shear 


We look now at the expression for the center of shear (7,,y,,/) as defined by Trefftz [10)}. 


The coordinates x, and y, are given by 


| . . 
1.=—F | feuds (3.1) 
7] | fords (3.2) 
I, ¢ al. 
FG 
Equation (3.2) may be rewritten as 
1 al ; aa 
Ys 61, | Jaxreda. (3.3) 
DF 


By Green’s formula we have 


— 
Ys 61, 


4 . , 
fp j2°nw— a? ds 
on 
( 


l , 2 . \ 
=-— fp [3a*n.e— a? (yn,—arn,)| ds. (3.4) 
61, . ; 
4 
The divergence theorem then yields 
| , | E ' 7 
Ust5y yd A 2] q rn gs. (3.5) 
—““UV /e “7 « 
DF a 
We introduce the notation 
: tf. i ai 
p Uetsy [awl (3.6) 
v “G 


and two harmonic functions h; and h,. Let 


6.= Phiznals, (3.7) 
C 
b.= 1p Oe ds= Phalyn.—zn,lds, (3.8) 
‘ on. P ‘ 
Cc Cc 
2 Ohe a 
bs Pi = Is (3.9) 
A 


bs 


J 


Combining (3.5) to (3.9) we obtain 





6-5, 4,+4. Di h,)(22n;—he)ds. (3.10) 
C 
Then by Schwarz’s inequality 
[b—b,—b.—b,|?< fre h, )°ds p( rn om ) ds. (3.11) 
si pedis. é - on 
c rf 
Again it is desirable to choose (g—/,) in such a way that 
PD (e—h,)ds=0. (3.12) 
C 
Using (2.23) and (2.39), we obtain finally 
L? £ 0 Oh, \ 4 Oh. \? 
[p—h, b.+5.)?< p. Tete as\ ri. =“) ds, 3.13 
[p>—bh,—b,+-b, Sa (3° a) XP 20, 9 (3.13) 
( ( 


We insert the known boundary value for gon the right-hand side and use the Rayleigh-Ritz 
. ¥ : A r 
technique to obtain close upper and lower bounds for y,. Note that if the boundary value 


problem 


Ah. 0 In Y 


oh, 


> rn, on (3.14) 
n 


can be solved, it is then possible to determine y, explicity without solving either the torsion or 
the flexure problem. 
In a similar way it is possible to obtain upper and lower bounds for the coordinate z,. 
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Half-Round Inductive Obstacles in Rectangular 
Waveguide 


D. M. Kerns 


(January 30, 1960 


Formulas are derived for the accurate calculation of the Jowest-mode, lumped-eiement 
representation of perfeetly conducting half-round inductive obstacles in rectangular wave- 
guide. These obstacles consist of either one or two opposed semicircular cylindrical indent a- 
tions extending across the narrow sides of the waveguide. They seem especially suitable for 
use as precise calculable standards of reflection or impedance in waveguide. Schwinger’s 


integral equation approach [1,2]! is used to obtain stationary expressions for the desired 
parameters as functionals of the surface currents on the obstacles. Upper bounds are ob- 
tained for one of the two parameters. Explicit formulas are derived for the values of the parameters 
under the assumption of n-term Fourier sine-series expansions for the obstacle currents. Rapid 
convergence is indicated by numerical evaluations for n=1, 2, and 3. In the process of 
obtaining expressions suitable for numerical calculation, an expansion (believed to be new) 
of the Green’s furetion of the problem is obtained and the sums of certain infinite series of 
Bessel’s functions occurring in this expansion are expressed in terms of definite integrals. 
\ brief numerical table of these sums, sufficient for the evaluation of the n= 1 approximation, 
is included, 


1. Introduction 


In this paper formulas are derived for the accurate calculation of lowest-mode, lumped- 
clement parameters for what may be identified as “single half-round” and ‘double half-round”’ 
inductive obstacles in rectangular waveguide. As shown in figure 1, the obstacles consist of 
semicircular cylindrical indentations extending across the narrow sides of the waveguides. 

The particular geometry considered seems especially suitable for obstacles to be used as 
calculable standards of reflection or impedance in waveguide. The geometry is well suited to 
electroforming, so that obstacles may be fabricated by this process as well as by machining. 
For obstacles producing standing-wave ratios of moderate values, the obstacle radius is large 
compared to high-standard machining tolerances, so that unduly close tolerances are not 
required, A change in waveguide cross section before and after the obstacle is avoided, elim- 
inating the need for waveguide components in odd sizes and permitting maximum flexibility 
of interconnection. ; 

In the present work the obstacle and waveguide surfaces are assumed perfectly conducting. 

The Waveguide Handbook [3] contains approximate formulas for semi-elliptical obstacles,’ 
which can be specialized to apply to the present problem. However, for the contemplated use 
in standards work, formulas with very much greater accuracy (10° or 10° times greater) are 
wanted, 

Brief descriptions of the present and some related work have been given previously [4}. 


2. Formulation of Problem 


We employ the customary complex electric and magnetic field vectors E= E(r), H= Hr), 
which satisfy Maxwell’s equations in the form 

V> E ' jovdd, 

VX H=jweE+ J, 


(2.1) 


Figures in brackets indicate the literature references at the end of this paper 
Che pertinent formulas are afflicted with a number of misprints; correct forms of these results for half-round obstacles are given later in this 


Pp per. 
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under the assumption of time-dependence represented by the (omitted) factor exp (jet). Here 
tis the time, w/(2r) is the frequeney, 7 is the imaginary unit, J is the complex current-density 
vector, and w and ¢ are respectively the permeability and the permittivity of the homogeneous, 
isotropic, nondissipative medium in the waveguide. The MKS system of units is employed. 
As shown in figure 1, we choose a rectangular coordinate system Oryz such that the in- 


terior of the waveguide is the space 0 r<a, O<y<b, o <>< o and the surfaces of the 
obstacles are given by the loci . 
I? =z + 2 0<r< R, 
(2.2) 
; 2 


Pe (a- r)?+- 27, 0<a—r< hy 


in the double half-round case and by 























Ré=77+- 2’, O<r<R (2.3) 
in the single half-round case. 
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Fiat RE lI. Single and double half-round inductive obstacles 


in rectangular wavequide 


The fields incident on the obstacles are, by hy pothesis, to be in the TE, mode, Traveling- 


wave forms for this mode may be written 


(2 4) 
H: -ysih (wre dae, + pr (wed ' cos (mraje le al 


E*+—sin (rr/aje*#e,, \ 

J 

The upper and lower signs refer respectively to waves going in the positive and the negative 
2-directions; the phase constant 8 is equal to [4°—(ra)]'?, where k= w(ue)!?; the wave-ad- 
mittance »=8/ (wu); and e,, e,, e, are the unit vectors of the system Ory2. 


r<hka<2r, so that the TE, mode, but no higher mode of the type TE, , will have real propaga- 


It is assumed that 


tion. 


We define “voltage” and “current” +,,(<,,),7, for the TE,, mode by means of the 


equations 
. 


st Dat Zan} sin (mr a), 
m ee. (2.6) 
H..=| "tm (2m) SIM (wT a). J 


Here the indices m=1,2 refer respectively to the “left-hand” side (+, R) and the “right- 
hand” (2, >) of the obstacle considered, and £,,,, //,,, are components (identified by the 
subscripts) of the TE,-mode part of whatever total electromagnetic field may be present in 
the waveguide at z=<,. From (2.6) and (2.4) it follows that the waveguide characteristic 


impedance, defined as the value of ¢,,(2,,)/7,(2,) for a pure traveling wave incident on the side 


m of the obstacle, is equal to unity for m=1,2. 
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The desired lowest-mode, lumped-element deseription of an obstacle is contamed im its 
impedanee matrix Z, which characterizes the linear relations 


/ 


. Zit / Z sols, t= Zt, Zool (2.7 


imposed by the obstacle on the terminal variables ¢,,,7,,. The values of the Z,,,, depend upon 
the location of the terminal planes z= 7, and z= 72; in the present problem it turns out to be 
convenient to extrapolate these planes to 2, 72.—0. 

Now, the equality 7), —Z.. follows from the structural svmmetry (and the symmetrical 
disposition of the terminal planes) with respect to the plane :=0; Zj.—Z.; is assured by both 
symmetry and reciprocity; and, since losslessness has been assumed, Z is pure imaginary. 


Thus Z is of the form 


IN IX, 
IN At 


avnd there are only two independent parameters to be determined. 

We do not obtain formulas for Y,, and Y,, directly; but rather, to exploit: the symmetry 
of the problem more fully, we consider the impedances obtaining under modes of excitation in 
which the electromagnetic field is either svimmetric or antisvinmetric with respect to the plane 

0. Thus in the symmetric ease, 7(0)—7,(0), 7 (0)— 7,10), and the impedance may be 


written 
Ace=Vie(0)/r,,,(0), (2.9a 


where we have added a subscript distinguishing the “even” ease. In the antisvinmetric 
odd’) case, 7,00) r(O), 7,(0) »(0), and the unpedance may be written 


Joo = Vil Odo O), (2 Q))) 
again adding a distinguishing subseript. The relations 
2yN = Zeet Zou 2)X12= 4.2 —4oar (2.10) 


follow direetly from (2.7), (2.8), and (2.9). ([t may be noted that the symmetric and anti- 
svmmetrie field distributions correspond to the eigenvectors of the matrix (2.8) and the quanti- 
ties Z,, and Zo are the eigenvalues of this matrix. 

An equivalent network for the obstacles is of interest and is presented here for convenient 


reference, For a 7T-network, using the sign conventions shown tn figure 2, we find 


4,\= Le, 242= Fee — bem (2.11 


where Z, and Z, are respectively the series and the shunt elements of the symmetrical 7. It 
may be remarked that the obstacles considered are called “inductive” beeause for small radii 


the important element, Z., is a positive reactance 








i ‘ 
o MN ? oNA 0 
Vy iy Z2 ip V> 
Gucen + 0 
2,3 Ze 2, = + (Zee - Zoo) 
Fiaure 2. Equivalent network. 
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3. Analysis 
3.1. Reduction to Scalar Form 


Maxwell’s equations directly imply the general differential equation for E. 
vx<0& E=k°E— jon, (3.1 
Which must be specialized appropriately for the present problem. E is moreover subject lo 
the boundary condition that its tangential components vanish on perfectly conducting surfaces 
In the present problem both the exciting field and the waveguide-obstacle structure are 
independent of y in the range 0<y<b, and it follows that the whole field, secondary as well as 
primary, will be independent of y in the same range. Further, since the «- and >-comlponents 
of E must vanish in particular on the y=0 surface of the waveguide, these components must 
vanish for all values of y involved. Thus the electric field is of the form Eos, - e,; eq (5.1 


(with J=0) reduces to 
v*o+-k-o=—0 in », 


and the boundary condition becomes 


o—) On a 


where S is the cross section and Cis the boundary of the structure ina plane y— constant (fie. 3 
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Kigure 3 Relevant geometry and coordinates for the single 


and double half-round proble Wis 


The magnetic field corresponding to EB of the above form is 


_( Ob Oo rs 
H | Jwu '( €,- e ) 3.20) 
O: O: 
The surface current AH non Cis thus given by 
K — —(jon)-'n-v¢e 3.2b) 
evaluated on ©, where mis the outward normal unit veetor on ©. Thus A has only a y-com- 
ponent and may be written as A(r,z)e, 
Inasmuch as the whole problem can be dealt with in terms of @ and A, the problem is es- 
sentially a two-dimensional scalar one. In what follows it will be convenient to refer to @ 
simply as “the field” 
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3.2. Integral Equation Formulation 


We now introduce a Green’s function T=T(s, 2, 2’, 2’) such that T is the field that would 
be produced in the unperturbed waveguide at the point (x, 2) by a unit current filament par- 


allel to e, at the point (2’, 2’). As a function of (2, 2), T° satisfies the differential equation 
- bl - bl 4 ao , r , ‘ ‘ 

vr+k i jou 6(4—2") 6(z—2") (3.3) 

obtained from (3.1) with J=é(*—2’) 6(2—2’)e,, where 6 denotes the Dirac delta-function. 


r further satisfies the boundary condition [T=0 for z=0, a and represents outgoing waves 
for -—2’ >. An explicit form for T is [5] 


r(z, 2,2’ ,2') wu(Ba)~' sin (wr/a) sin (wr’ /a)e~* 
(jou a ) Say; ' sin (irr /a) sin (nmr’ /aje~ em FF, (3.4) 
where a, = y (ur a)?—F* for n>2 and the other symbols are as already defined. 
It is convenient to put 
r -wp(F'+-7G), (3.5) 


where / and G are real. We observe that F can be written in the form 


F'= (Ba)~'[y, (4, 2) ¥.(2’, 27) +-Wola, 2) Wola’, 2”), (3.6a) 
with the definitions 
¥,(4, 2)=cos Bz sin(me/a) 


; (3.6b) 
WY) (4, 2)=—sin Bz sin(rr/a). 


The functions y, and yp are respectively symmetric and antisymmetric with respect to re- 
flection in the plane z=0. 

With the aid of the Green’s function we can set up the fundamental equations of the prob- 
lem in integral-equation form. Let ~=y(s, 2) denote the exciting field—i.e., the field, due to 
remote sources, that would exist in the absence of obstacles. Then, with an obstacle present, 
the total field @ can be written as the sum of the exciting field and the secondary field due to 


the current on the obstacle: 


$(2,2)=¥(2,2)+ | T(2,2,2’,2’)A(2’,2’) dC’, (: 


ef 


~I 


~~ 


where AY is the surface current on the obstacle and the integral is a line integral going over the 
contour (, of the obstacle (Cy) consists of two parts in the case of the double half-rounds). If 
the point (r, 2) is taken on (Cy, @ must vanish and we have 

O=—Yr,2z)+ I'(z,2,2’ ,2’ )K(2’ ,z’) dC’, (7.2) on Ca. (3.8) 


e Cc 


This is an inhomogeneous integral equation of the first kind determining A for a given y. The 
equation reduces to slightly different forms in the symmetric and antisymmetric cases, which 
we now consider. 
a. Symmetric Case 

Let the exciting field be Ay,(2, 2), where A is an arbitrary amplitude and y,, as defined in 
(3.6), is a symmetric function of 2. The surface current on the obstacle will then also be a 
symmetric function of 2; we denote it by A,. We now examine (3.7) for 2< <2’, say, to find 
the lowest-mode component and thus to find ¢,,(2) and 7,.(2). For << <2’, (3.7) becomes 


o(2,2)—A cos Bz sin (rr/a)—(Ba)~'wyu sin (rr/a)e* | ¥.K, de, (3.9) 
Co 
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since integrals of the type | YX, dC vanish because of the mutual orthogonality of syvm- 
C, 


metric and antisymmetric functions. By comparison with (2.6) it is seen that 


v, (2) =A cos Bz—(Ba)~ wy oe} | VK, de. (3.10) 
Cy 


Using (3.2a) to obtain the z-component of the magnetic field associated with (3.9) and again 


referring to (2.6) one finds 


i, (2) =—jaAl sin Bz+ (Ba) ~'wp « | VK, dc. (3.1 1a) 
( ) 
Equations (3.10) and (3.11a) hold directly from the definitions for 2: Pe and by extrapolation 
for z=0. In particular (3.11a) vields 
t1,-(0) = (Ba)-' wn | VK, dC. (3.11b) 
Cy 


Next, in the integral eq (3.8), we separate the Green’s function into real and imaginary parts as 
in (3.5), use the orthogonality property of even and odd functions, and thus find 


W.(4,2) %y,-(0) =Jop Gizs.2.% 2) H(z’ ,2') dC". (3.12) 
Ca 
This equation and (3.11b) together furnish a definitive mathematical statement of the svm- 
metric part of our problem. 
b. Antisymmetric Case 


In this case we let the exciting field be Ayo(z, 2), Where A is again an arbitrary constant and 
Y, defined in (3.6), is an antisymmetric function of 2. The current, Ay, on the obstacle will 
then also be antisymmetric with respect to 2. In the same manner as in the symmetric case 
one finds, for z<0, 
t,0(2)=A sin Bz + (Ba)~"jape** Wh dC. (3. 14a) 
J Cy 
In particular, 
r, o(0)= (Ba) ‘jen | YK dC. (3.14b) 
J Co 
Further, 
i1.o(2)=jA cos Bz—(Ba)~"jope | WAAC. (3.15) 
. Cy 
and the integral equation becomes 
Wo(r,2)%1,0(0) Wu G(r,2,2" 2’ )Ko(2’ 2’ )dC’. (3.16) 


a 


This equation and (3.14b) together furnish a definitive statement of the antisymmetric cases 


of our problem. 


3.3. Summary of all Cases 


In the case of the double half-rounds the above integrations go over a two-part contour 
consisting of an upper semicircle (, and a lower semicircle C;. An equivalent problem involving 
only one of the parts of the contours, say C;, can be formulated using the fact that the exciting 
field as well as the structure is symmetric with respect to the plane z=a/2, so that the current 
on the obstacles must also be symmetric with respect to this plane. This symmetry is quite 
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independent of the symmetries with respect to 70 distinguished by the subscripts “e” and 
“9” and so these subseripts will be dropped for the moment. Thus, it is clear that 


|. VA dC | WK dC l. WK dC =2 P WK dc. (3.18) 


Further, 


| G(z,z,2',2’) K(x’ ,2’) dC’ | G (z,z,2’,2’) K(xr’,2’) dC’, (3.19) 
Cc. Jo 
where G™ is defined by 
G@ (z,2,2’ .2’)=G(2, 2,2’ ,2’)+Glz,z.a—r’.2’). 
Equations (3.18) and (3.19) enable the desired restatement of the double half-round problems. 


In what follows, superscripts 1,2 will be used when it is desired to distinguish quantities 
associated with the single and the double half-round problems, respectively (this means in 


particular G=G°). Integrals will be indicated by means of the convenient scalar product 
notation, e.g., 
(yp, Kk) | VK dC’, (3.20a) 
Cc) 
and by the operator notation 
Gh | G(z,2,2" ,2’) K(z’,2’) dC’. (3.20b) 
« C; 


As indicated, all such integrals are to be taken over (,, the lower semicirele. As a further 
notational convenience, we introduce the “normalized” surface currents, 


T,=jopk,/v,,(0), a wk, /7;..(0), (3.21) 


for both the single and the double problems. 
Our results thus far may now be summarized as follows. For the symmetric cases we have 


y= G I, 1/X.” =h(y,,1.)/(Ba) (3.22) 


and for the antisymmetric cases 
y= GJ, Nye = —h(Wo,Te)/(Ba), (3.23) 


where h=1 or 2 for the single or the double half-rounds, respectively, and X,,, Xoo are the 
reactances corresponding to the (pure imaginary) impedances Z,,., Zo. defined in (2.9). 


4. Solution for the Reactances 


4.1. Application of Rayleigh-Ritz Method 


In each of the problems specified by (3.22) and (3.23), an unknown reactance or suscept- 
ance is proportional to a scalar product 


M=(y, D, (4.1a) 
where the function J is determined by an integral equation, 
GI=y, (4.1b) 


in which y and G are given and G@ is symmetric. From these equations one may easily con- 
struct the “stationary representation” [1, 2] 
M=(y¥, D?/U, GD) (4.2) 
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for M. This constitutes the starting point for the method of Rayleigh and Ritz [6], which 
we wish to use to obtain an approximate expression for 17. One assumes in (4.2) an approxi- 


mation for the surface current in the form 
l.=Z 1, (4.3) 


where the 2, are coefficients to be determined, the /,; are members of a suitable set of basis 
functions (to be chosen explicitly), and summation from 1 to n over repeated indices 7, /, a 
is understood. Expression (4.2) becomes an ordinary function of the x, and is to be subjected 
to the conditions of stationarity 0.M/dr,;—0, i—1, 2, ..., n. This leads to the system of 


equations 


Ga (w.¢,/.M, )e 


where 
¢;:=(¥, AN G (fi, Gi , (4.4) 


and \V, denotes the now-determinate approximate value for V/. If the arbitrary normalization 


of the x, is chosen so that x,c,—.M,, then 
e (G 1) js 


where (G~'),, denotes an element of the inverse of the 1 <n matrix of the G,,, and it immediately 


follows that 


M,=(G-"') se €;. 


Finally, we write this quadratic form as a ratio of determinants, 


0 C) i 4 C, 
C1 Gy a ae 2 Gi, 
Cc, Gr, - — Gr, , 
—M, (4.5 
G); a a a (i, 
G l . =e ie Gi, ’ 


Thus JV, is expressed in terms of the known quantities (4.4); if the objective is only to calculate 


the reactances, it is not necessary to caleulate the 2 
4.2. Definiteness of G; Upper Bounds for Xv 


Although in our problem G is not a definite kernel, it does have definite character with 
respect to a suitably restricted class of admissible functions. In fact, @ is positive definite 
with respect to functions that are antisvmmetric with respect to 2 and vanish for 2 >A,/4, 
where \,=27/8 is the “guide wavelength” of the single propagating mode.’ In terms of the 
integrals of interest here, this means that if T, is any assumed surface current distribution, 
continuous, not identically equal to zero, and antisymmetric with respect to z, then 


‘a GT,) -(). (4.6) 
provided R<— ,/4. 

Now, for a positive definite kernel it is well known that the right-hand side of (4.2) is not 
only stationary, but also a maximum for the true 7. For this result we offer the following 
brief proof based directly on the Schwarz inequality.’ Assuming @ to be positive definite, 
we wish to compare the approximate value for M/ given by 


W 7. y)° CT. GI), 


3 This result is suggested by qualitative physical considerations and has been verified analytically by the author. It has also been verified 
that @ is not definite with respect to functions symmetric with respect to 

4 Proofs of this result in the literature known to the author, e.g., [1,2], employ the so-called bilinear expansion fer the kernel involved and do 
not vield the “‘only if’’ part of the result (4.7). 
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where / is an admissible approximation for the surface current, with the true value given by 
(4.1a) (as well as by (4.2)) when J satisfies (4.1b). From (4.1b) it follows that (/,~)=(/, GJ); 
the Schwarz inequality for positive-definite symmetric transformations [7] gives 


(7, GI)’ < (7, GT) (1, GI). 
Since (J, G1) =(1.W) =M, we have the result 
WeM. (4.7) 


Moreover the sign of equality holds in the Schwarz inequality and hence in this result if and 
only if 1 is proportional to J. 

Clearly, (4.7) is applicable when G@ is merely positive definite with respect to all functions 
that need be admitted; by (4.6) this will be the case in the antisymmetric parts of our problem 
provided R—),/4. Hence, subject to this inequality, in the antisymmetric cases we shall have 
M,<M (we do not anticipate an exact result for finite n); this in turn implies that numerical 
values for Xo will be (algebraic) upper bounds for the true values—assuming, of course, that 
numerical evaluations involved are sufficiently precise. 

Unfortunately one can not obtain a bound for X,, in this simple manner. Furthermore, 
derived quantities of practical interest (such as (4.17)) will usually depend upon both X,, and 
Now. Thus in general no bounds for such derived quantities are determined. 


4.3. Basis Functions 
We introduce plane polar coordinates 7,@ such that 


r=r sin @, “=r cos 6. (4.8) 


The path of integration (, is then given parametrically by 
rR sin 86, -=R cos 8, (4.9) 


where R is the obstacle radius and 0Zé@<Zzr. On this path the surface current J=J(R sin 8@, 
R cos @) becomes a function of 6, which we denote simply by /(@). We observe that sym- 
metries with respect to 20 are equivalent to symmetries with respect to 6=7/2. 


As basis functions we take 
f;:=2(rR)™ sin (27—1)0 (4.10) 
for the symmetric cases (expansion of /,), and 
f,=2(xR)™ sin 210 (4.11) 


for the antisymmetric cases (expansion of J)). The two sets of functions are complete for the 
expansion of symmetric and antisymmetric functions, respectively, in the interval (0,7). The 
completeness gives good assurance that for sufficiently large n, ., will approach arbitrarily 
close to the true VM. For practical numerical calculation it is vital that the convergence of M, 
be rapid. In the present instance rapid convergence might be anticipated on the grounds that 
the true J, and J, may be expected to be smooth (possibly infinitely smooth) functions of 6 
vanishing at 0 and x, so that their Fourier sine series should converge rapidly.* 


—_—_—_—_— 


' The points (R, 0) and (R, ) are the vertices of the corners formed where the halround centered at O meets the waveguide wall. The van 
ishing of the surface current in these vertices is equivalent, by (3.2), to the vanishing of the components of the surface magnetic field normal to 
the vertices. That these field components should vanish in the limit as a vertex is approached follows from properties of Sommerfeld’s solution 


114] of the proble m of diffraction by a wedge of arbitr ivy angk 
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4.4. Evaluation of Integrals 


Integrals of the type (y, /;) may be obtained relatively easily with the aid of the well-known 
expansion [8] 


it sin ~ ‘ jne 
e” wa > J,.(¢)e*, 


i= 


wherein J, is the n™ order Bessel function of the first kind. In this expansion we put ¢=kR 
and a=+60—x+(/2), with k sin x=z/a, k cos x=8, and employ (4.9). The expansion 


e sin (rr/a)=2>5 59" 1S, (kR) sin n9 sin nx (4.12) 


is then easily obtained. The integrations vielding the desired quantities are now elementary; 
one obtains 


rp - 2(—)” 1) 2 J,(kR) sin PX; pP 21- w (4. 13a) 
for the symmetric cases, and 
¢,=2(—)?-?? J,(kR) sin px, p=, (4.13b) 


for the antisymmetric cases. (It is encouraging that these quantities decrease rapidly with 
increasing index. ) 

The evaluation of the integrals of the type (/,, Gf;)=G,, is considerably more complicated 
than the evaluation of the ¢;, and the task has been relegated to the appendix The results 
obtained may be summarized conveniently in the following form, which comprehends four 


cases: 
Gi" ] 
to [p,(AR)5,, +0," (ka)] (4.14) 


Here s,=sin px; p,(kR)=Y,(kR)/J,(kR); in the symmetric cases, p=2i—1 and g=2j;—1; in 
the antisymmetric cases, p=2i and g=2); the functions o,") are defined in the appendix (and 
tabulated briefly in table 2); }°, denotes the Bessel function of the second kind; and all other 
symbols involved have been defined previously. 


4.5. Results 


The value of a reactance element corresponding to an n-term expansion (4.3) will be 
called an “n™ approximation ‘ 


and will be written simply?) or YXoj’ without special notation 
indicating the value of n being considered. An explicit expression, in determinantal form, for 
X,? in the n™ approximation may be obtained by combining (3.22), (4.5), and (4.14); viz, 


Piton O\ on-1 
h) h - 
yo fe a; Qn 1 ~~ <4 Peon oe Ton 1,2n 1 (4.15) 
ae 
4h \0 s, Son—1 
S + h 
I Piro Oi on-1 
h) 
Son 1 F\ 2n 1 P2n l Fon 1,2n-—1 
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(In this and in the following equation the arguments kR and ka are understood.) Similarly, 
(3.23), (4.5), and (4.14) yield 


0 &, Sen 
(h (h) 
2 Pot 9’ F2.9n 
> (h) (h) h) 
yn __ 4h Son 2. 2n e Pon + Fon, 2n\. (4.16) 
«oo . y th th) 
Ba Pot Oo9 “a To an 
h) (h) 
G9, 2n */-* 2 Pe2n T Ton, Qn 


The determinantal form of these expressions seems convenient enough for practical calcula- 
tions with » at least as large as 3. However, it may be remarked that, because of the special 
form of the diagonal elements, the determinants can be expanded in the same way as certain 
determinants occurring in the Fredholm theory of integral equations. This expansion may 
be convenient for the consideration of arbitrarily large n, but is not needed here. 

Once the reactances Y,, and Xo are determined, other quantities or parameters associated 
with the representation of an obstacle are of course also determined. In particular, the “volt- 
age’’ standing-wave ratio seen on one side of an obstacle when the other side is terminated 
in a matched (reflectionless) load is a familiar quantity having immediate physical significance, 
and it is convenient to discuss some of the results in terms of this quantity. Recalling that 
waveguide characteristic impedances were chosen equal to unity and referring to the equivalent 
network, figure 2, it is easily found that the standing-wave ratio in question is given by 


n= (1+ |S|)/(1—[S]), (4.172) 


where |S), the magnitude of the associated reflection coefficient, is 


4 r  \ 2-1/2 

[S| [ (er) | (4.17b) 

In order to get some idea of the behavior of the sequence of approximations, the reactance 

elements and the related VSWR (4.17) have been calculated for n=1, 2, and 3 and for a range 

of values of ka and kR. Table 1 presents more or less typical results in terms of the values 

of VSWR. The convergence of the sequence of approximations, judged on the basis of 

numerical results exemplified in the table, appears to be very rapid; by the same token, the 
first approximation furnishes a rather good result for VSWR’s up to about 2. 

Thus far it has not proved feasible to determine the manner of convergence of the sequence 

of approximations by theoretical means. In view of the character of the particular problem 

at hand, it seems probable that the numerical evidence may be relied upon. 


TABLE 1. Convergence of VSWR 


ka=4.5 
Approxi- VSWR, VSWR, 
kR mation Single half- | Double half- 
round round 
| 1 1. 0370968 1. 0776492 
0.2 2 1. 0370970 1. 0776499 
| 3 1. 0370970 1. 0776499 
| 1 1. 4547668 2. 8411581 
0.7 2 1. 4554655 2. 8416266 
| 3 1. 4554655 2. 8416268 
| 1 2. 1070625 15. 900756 
1.0 2 2. 1125072 16. 009131 
| 3 2. 1125112 16. 009479 





Note: Maximum computational inaccuracy in 
the above figures is estimated to be approximately 
+1 in the sixth significant figure; eight figures, as 
given, may be significant in indicating the behavior 
of successive approximations, 


123 








vSwR 


Further calculations have been made regarding ka and R/a as independent variables. 
This choice of variables corresponds to the practical situation where one has a given obstacle 
(having a given value of R/a) that is to be operated at various frequencies. Figures 4 and 5 
show the matched-termination VSWR (4.17) resulting for various values of R/a and Shes 
where / is the operating frequency and f, is the cut-off frequency for the TE, mode (and 


Sife=ka/r). (Precise tables (based on the third approximations) giving basic parameters and 










































































possibly some derived quantities are being prepared for separate publication. ) 
mn . . ° > . . 
lhe first approximations appear to be of useful accuracy for some purposes and represent 
an appreciable extension of previously available results. Consequently, a brief tabulation 
of the o-functions, sufficient for the calculation of first approximations, is presented as table 2. 
‘ ° ¥ ° ° - . , . 
For convenience of reference we write down (4.15) and (4.16) for the special case n= 1, 
! 
_ rT. |) lixe pa > h 
| | | | | X ce wags? 3 lpi(kh) tai (ka)], (4.18) 
t + + } t } + 4 { Shs? 
R/a=0.28 
3 O-—_++ + + + + + + + + + + — 4} 2 
iNo5 
} > (h) 2 > h ‘ 
0.24 2 va A 00 [po(kR) tT Fo9 (ka)| 2 (4.19) 
~<a OS Sb ee ed on oo Ba 
0.20 
}—___+ - + + + + + + + . ao = . 2 
Pos | | and note that s,;=x/(ka) and s,=278/(k’a). 
a ee | 4 } } 4 } } 4 4 ecusaiel 
Sa eS SS eS eee ee ee ee TABLE 2. Values of o%)(x) 
ate a 
20 + + + + + + + + + —— . 
0.16] | oh) a, 1 5 (r 
-——_+ + + + + + + + + + —~ | — - a 
| 3.80 | 2. 055435 0. 3170852 3. 960997 ~2 544410 
= ame ws 7 T T T T T T T a. 3. 85 1. 976710 1982432 3. 857142 2 
2 3.90 1. 899606 —. (90TS960 3. 756554 ~1.4 
— S + + } } } Ps } } } } | 3.95 1. $23905 006197600 3. 659060 —1. 633535 
9.08 | 4.00 1. 749764 09356160 3. HASLO — 1. 376168 
5 + + 4 4 + Ps 4 —_ j 4.05 1. #76799 1720644 3 5 1. 138910 
04 — 4.10 1. 604988 2424036 3. 383659 —(). 9201812 
10 — + i 1 | 4.15 1. 584222 3052036 3. 297087 —. TISA5R0 
1.2 . r + 
‘ 4.0 1. 464400 SH1L0476 3. 212923 5327088 
4. 25 1. 305414 $104.54 4. 131048 3614372 
. rou ue f . 4. 30 1. 327159 PSOE 3. 051353 2036444 
k IGURE 4. VsSu R versus J J- jor single half-rounds. 4.35 1. 250533 . 4918312 2. 973735 —. 05831720 
4.40 1. 192430 5246272 2. 8OS004 O7546600 
4.45 1. 125746 5526568 2. 824341 . 1985600 
4.50 | 1. 050372 STH2488 2. 752387 3117472 
4. 5S 0. 9931968 5USHO4 2. 682146 4157360 
4.60 9271072 HL 12796 2. 613542 5112028 
4.65 SHOUSSH 6232416 2. 46500 5OS7 500 
4,70 . TOAGORS 6318024 2. 480046 OTRUBIS 
4.75 7281208 6371680 2 416812 . 7522756 
4.80 HOLLO76 6395192 2. 354034 S192504 
4.85 SYS 5004 6390196 2. 292551 . SSOBTSO 
4.90 251500 6358 10S 2. 232301 Y35R8506 
4.95 4558572 6300236 2. 173228 YS62656 
5, 00 BN4276 H217692 2. 115279 1, 031884 
5.05 3136376 6111440 » O5S400 1. O7 3088 
5.10 2402376 SYNZIS4 2. 002545 1. 110027 
5.15 . 1649416 SS30R36 1. O47661 1. 143125 
5. 20 08742720 . WSTH2Q8 1. 893705 1. 172504 
§. 2% . 007320000 H2044 1. 840631 1. 198667 
5. 30 O75S81400 5246072 1. TRABOT 1. 221571 
5 35 1624856 SOONG 1. 736963 1. 241508 
5. 40 —. 2533004 - 4750804 1. 686288 1. 258672 
5.45 | S4SUSS4 HH9012 1. 636332 1. 273236 
5. 50 1504328 . 4168292 1. SS7059 1, 285370 
5. 55 —. 458722 . BRBSY12 1. 538434 1. 295225 
5. 60 ~. 6752200 3486376 1. 490420 1. 302940 
5.65 SOT6548 . 3106864 1. 4429081 1. 308654 
5. 70 —. 9402032 2697056 1. 306088 1. 312488 
§.75 —1. 094004 2OSH4ASS 1. 349701 1, 314555 
5. 80 1. 266887 1778340 1. 303798 1. 314964 
. - — - . Note: Accuracy of the above figures is estimated to be ipproximately one 
Fieure 5. VSWR versus f/f. for double half-rounds. | part per million, . 
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For sufficiently small #2, the expressions just given may be simplified still further. For 
small &P, the quantities p,(AR) and p(kR) become large, in accordance with the formula 


(ER) ce — ot (e—1)1(-2-) 51) 
Pa( ‘ ~— | = ° ( kR ) (q> ’ 


and thus determine the qualitative behavior of Y,, and Yo. Using this approximation for 
the p’s, neglecting the o’s, and introducing the guide wavelength \,=2z2/8, one obtains 


am. mfey ; 
X. in, (cr) ' (4.20) 
4 
xX” — (=). (4.21) 
g 


(These are the correct forms of the results in the order of approximation of the more general 
formulas in the Waveguide Handbook |3).) Although these formulas give a rather good picture 
of the behavior of the reactances, they fail to reflect the interesting qualitative feature that 
Xx“? and VY“) must eventually become equal to each other as 2 approaches a/h (as can be seen 
physically). 


5. Appendix 


5.1. Expansion of Green's Function in Terms of Cylindrical Wave Functions 


In order to arrive conveniently at the respective types of Green’s functions needed for 
the single and the double half-round problems, we consider a waveguide of width w (instead 
of a), bounded by the lines s=0 and s=w; and, in addition to the electric boundary condition 
at s=—0, we consider both electric and magnetic boundary conditions at r=w (by “electric’’ 
or “magnetic”? boundary conditions it is meant that the Green’s function or its normal deriva- 
tive, respectively, vanishes on the boundary). The radius vectors of the source-point and the 
field-point will be denoted by r’ and r, respectively. 

Now, it is well known that the two-dimensional “free-space’’ Green’s function, satisfying 
(3.3) and representing outgoing waves at infinity, is —(wu/4)//)(kir’—r}|), where 77, denotes 
the Hankel function of the second kind (as is appropriate for exp (jot) time dependence). 
It is also well known that the desired Green’s function, ®, satisfying (3.5) and the required 
boundary conditions at s=0, r=w, and at z= +, can be obtained in the form of a sum of 
free-space Green's functions by the method of images [9]. In fact 


> (wu/4) >) s"[(F(kipi—r|)—Hy(kip;—r))], (5.1) 
where p* = (2nw+-2’ )e,+-2’e, and s=+1 or —1 according to whether the boundary condition 


at sw is electric or magnetic. By inspection of the array of sources (fig. 6) it can be seen 
that this function has the proper symmetry to satisfy the boundary conditions: it is anti- 
symmetric with respect to s=0 and it is antisymmetric or symmetric with respect to z=w 
according to whether s= +1 or —1. 

Equation (5.1) can be transformed into a sum of elementary wave-functions of polar 
coordinates by repeated application of Graf’s addition theorem for Bessel functions [10]. 
We introduce polar coordinates (7, 4), (7”, 0’), and (pp, @,) for r,r’, and p, in the manner shown in 
(4.8), and we assume r and 7’ both to be less than w, so that p, >w for n # 0. Considering 
the triangles with sides p,, r, and p,—r and applying the addition theorem, we obtain 


—4(wu) ‘b= 9 Jy (ker)H,(kr’) 
pon Ip(kr’ )H, (kr) 


¢ 1P9(¢ jper eiPe") 


+> 8" D3 J, (kr)e”"(H,(kpte-"t —H,(kpa)e-"*a], 
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Figure 6. Source points in the method of images. 


Relative signs of sources are indicated by the adjacent symbols +, 
+s, and —s; the role of s is defined in connection with (5.1). 


where the terms arising from n=0 appear in the first line and are to be omitted from the primed 
summation in the second line; the upper and lower alternatives in the first line apply according 
as rr’ or r>r’, respectively. Again, by the addition theorem, 


x 


H,(kp*)e~ n= >> « g(kr’ )H,4,(\2nkw)) exp [7(x/2) (q—p) sgn n+) 90’): 
ef 
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Combining the last two equations and rearranging terms, one finds that (5.1) can be written 
in the form 
a ae , J (kr) H,(kr’ ) 

p wu >. >* sin pé sin gé’ ghee 5, 

24 2 sin} YJ, (kr’)H, (kr) °?8 

+2J,(kr)J,(kr’) cos |(x/2)(p—q)] >) s"|H,-,(2nkw)— H,.<(2nkw)| }- (5.2) 
n=1 

If now one puts w=a and s=+1, this expression becomes the desired representation of T of 
the text; similarly, if w=a/2 and s=—1, (5.2) becomes the desired representation of a Green’s 
function T™ related to G® in the same way as T is related to @ (ef.3.5). Thus, employing the 
relation (3.5), the relation 7,=J,—jY,, the notation (3.20), the basis functions (4.10, 4.11), 
and carrying out the integrations with respect to 6, 6’, one obtains 


9 \2 
(~,) (sin pé, G@™ sin g0’)=—J,(kR)Y (kR)b,,—T (RR) FT (RR) cos [(x/2) (p—Q of (ka), (5.3) 


where the functions of) are defined as 


oy. (ka) =2 24 [Y ,-,(2nka)—Y ,4,(2nka)], (5.4a) 


a}? (ka) =2 >) (—)"[¥ p-,(nka)—Y ,,,(nka)). (5.4b) 
n=1 
We note that the above expressions are needed only for positive integer indices p, q that 
are both even or both odd. We observe also that the form of the expressions is such that 
those with unequal subscript indices can be evaluated simply in terms of those with equal 
subscript indices. Further, in view of the easily-verified relation 


a}? (ka) =20})) (ka) — oS) (ka/2), (5.5) 
values or formulas for the o’s “‘of the second kind” can be derived readily from the same for 


, oe . . * rr . yth) . . . . 
the o’s ‘of the first kind. Thus, the evaluation of the Gy, which, on the face of it, originally 
involved double integrals with ka, kR, p, g, and h as independent parameters, has been reduced 
essentially to the evaluation of sums of the form 3) Y2,(mzr) with m<a<2m or, 24r<2<4m and 

n=l 
s=0, 1, 2,.... It should be observed, however, that these sums converge extremely 
slowly: for the accuracy desired in the present circumstances the number of terms required 
would be of the order of 10° or more, at least for the most important o’s. Needless to say, 
some method of evaluation other than direct summation is indicated. 


5.2. Evaluation of the Bessel-Function Series 
Two methods of evaluating the series defined in (5.4) will be discussed briefly here. 
a. Conversion to Definite Integrals 


An evaluation of the desired series may be obtained in terms of definite integrals as follows. 
We start with the series of Hankel functions (of the second kind) 


S)(2) > H(nz), (5.6) 


n=1 


where « is real, positive, and not an integer multiple of 27. For the Hankel functions we 
use the integral representation [11] 


1 . 
Hy(x)= | en ir sing+ Me 
Tv + 
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where ¢=£+-jn and C is a suitable path, to be chosen in a moment, between the limits —je 
andaz+jo. With the aid of the well-known expression for the sum of a finite number of terms 
of a geometric series, one finds that the \/-term partial sum of (5.6) can be written 


M . IM] . jMr sin ¢ 

S07] .) l _¢ dé Pe ( Me. sa 
Aid jr sing jr sing ¢ ag (9.4) 

n=1 TICE — ‘ l Tce : . ] 


Among the infinitely many poles of the integrands in these integrals, we note in particular 
e land t —r—t., where &, is the principal value of sin ~'(2mm/r), m=0,1,2, . . ., GY, and Q 
is the largest integer such that 2Q2/r<1. The path C is now chosen to comprise (except for 
indentations) the negative imaginary axis, the real axis between 0 and z, and the upper half of 
the line r+-jn; indentations, of radius p, are made at each of the above identified poles in the 
manner shown in figure 7, which illustrates C for the case Q=1. Taken on this path, the 
second integral in (5.7) approaches 0 as .V/ increases indefinitely; as a detailed examination 
shows, this happens because on the various parts of the path the integrand approaches 0 or it 
becomes increasingly oscillatory (or both) as 1-0. Hence we have 


1 r Md er 
S\ (r)=- | oa 
‘ ~ 


. 
Tr jrsing _ l 


This result can be expressed in terms of real integrations plus contributions from the 
indentations. In the limit as p-—>0, the indentations vield 


, ad é - ) 
PY. J R* x3 (5.Sa) 
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JXE jrx(w—t 
E78 m ” ¢ Sm a 
# , om “ (5.8b) 
y r°— (2mm)- y 0° —(2mn)- 


Here the notation corresponds to that used for the respective poles and the factor 1/ standing 
in front of the integral is included. Provided \ is an even integer, the straight-line segments of 
the path of integration yield 


Be ee } "= cosh An d ee | 1 ; i 
a 6.,+lim J l2 | ee sinh a *—P | COSPAE cot (5 v sing ) dé |, (5.9) 
‘ /?P ‘ af - ‘ 


“ p0 T, 


where 6,, is the Kronecker delta. It should be noted that when the integrand in the finite 
integral is singular, the Cauchy principal value is to be taken, as indicated by the symbol P 
(the indicated limit as p->0 essentially is also a Cauchy value). 

Upon comparing (5.4) and (5.6) it is seen that 


a}, (ka) =2 Im [S,, (2ka)—S, (2ka)|. 
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It happens that for \=2p (where p=0, 1, 2, . . .) the net imaginary contribution of (5.8) to 
S,, vanishes. Consequently (5.9) vields 


°F 


‘sin? pé cot (ka sin &) dé (5.10) 


2ka sint 
T ” 


' 8 ( * sinh? pndy +4 P| 
J 0 T 


J 0 


wherein the integrands are now well-behaved at their lower limits and thus we were permitted 


to put p=0. The corresponding result for o,, is 
i 8S (°° sinh? pndyn 4, (7? ., ka. rs 
o}?) (ka) —=|. fa stake t—=Pi sin? pé tan ( 5 sin Ejdé, (5.11) 


which is easily obtained from (5.10) with the aid of (5.5). 
(An interesting by-product of the above development, obtained by taking the real parts of 
(5.6), (5.8), and (5.9), is 
cos [2p sin”! (2 mm/z)} 
, 


. l @ 
Soot DoJ2p (nx) r +2>>- 


m=1 “f (2m)? 


tol— 


where p=0,1,2,. . .; cand Q are as specified above. This result was obtained by Ignatowsky 
(12) by a different method.) 

Thus far numerical results for the o’s have been obtained primarily from (5.10). The 
integrals were evaluated numerically by Paul F. Wacker and William W. Longley, Jr., using 
Legendre-Gaussian quadrature formulas programed for machine computation.® Values 
for of) were calculated from values of o{) using (5.5). Results for the two lowest values of 
p are given in table 2. (More extensive tables are planned for the previously-mentioned 


separate publication. ) 


b. Conversion to More Rapidly Convergent Series 


An alternative method of evaluating the o’s is afforded by a transformation of the Bessel- 
function series involved into more rapidly convergent series. This was accomplished a good 
many years ago by W. von Ignatowsky [12] in connection with his work on diffraction by 
gratings—a subject that will be recognized as being rather closely related to the subject of the 


present paper. Ignatowsky obtained 


: 9 Q | - | 
> 2K K a 
«> Y,(nz)=log “4-33 -+ 3 ( en ee: (5.12) 
n=1 Y n=1 n=Q+1 NM a 1 . 
—_— 1 12 (2«)*(p- IIb, 
us > & ) op( Nr) ») ao pe » f = ' = 
n=\ “ 2 p=1 (2u)!(p bh): 
@. sin |2p sin” '(nx)| (nxk—y n*«*— 1)” 
2 ~ ? x ~ 
—o Be a (— Pc 2 (5.12b) 
n=] y |—ner n=Q+1 yun —l1 
In these expressions k=2a/r; p=1,2,3, . . . ; Q@ is as defined above; log y=0.57721566 . . . is 


the Euler-Mascheroni constant; and the B,, are the Bernoulli numbers, here so labeled that 
B,=1/6, By=1/30, By=1/42, Be=1/30, By =5/66, ete. Equation (5.12a) is fairly well known 


at the present time. (A few numerical values are given for the series >) (—)"Yo(nx) in 


n=l 


6 The estimates of accuracy of computed quantities given in this paper were made by Mr. Longley, who was responsible for most of the detail 


of the calculations 
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[13].) The series for o{? obtained from (5.12a) and (5.12b) has been duplicated by an in- 
dependent mede of derivation; this represents a partial check of the complicated expression 
(5.13b). Thus far in the present work, Ignatowsky’s series have been used only to spot check 
some of the numerical results obtained from the definite integrals. The series are especially 
suitable for either estimation or calculation of the o’s when p is large. j 


The author has benefited from the discussion of many points with Paul F. Wacker and is 


indebted to David F. Wait for assistance with some of the theory and computations. 
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Selected Abstracts 


4 comparison of experimental and theoretical rela- 
tions between Young’s modulus and the flexural 
and longitudinal resonance frequencies of uniform 
bars, S. Spinner, T. W. Reichard, and W. E. Tefft, 
J. Research NBS GAA, No. 2, 147 (1960). 


The relations between Young’s modulus and the mechanical 
resonance flexural and longitudinal frequencies of two sets of 
steel bars has been established by an empirical method. 
The agreement between theoretical and experimental curves 
for flexural vibrations appeared somewhat better than that 
for longitudinal vibrations. 


On the mode theory of very-low-frequency propaga- 
tion in the presence of a transverse magnetic field, 
D.D. Crombie, J. Research NBS 64D, No. 3, 265 
(1960). 


The effect of a purely transverse horizontal magnetic field on 
the propagation of VLF waves is considered. It is shown 
that the magnetic field introduces non-reciprocity, and that 
for propagation along the magnetic equator, the rate of 
attenuation is less for west-to-east propagation than for east- 
to-west propagation. 


On the theory of reflection of low- and very-low- 
radiofrequency waves from the ionosphere, J. R. 
Johler and L. C. Walters, J. Research NBS 64D, No. 
3, 269 (1.960). 


The rigorous application of the magneto-ionie theory to the 
calculation of reflection coefficients for a sharply bounded 
ionosphere model is discussed. This is a comprehensive 
paper on the details which pertain to the rigorous evaluation 
of the reflection coefficients. The paper is illustrated with 
computations applicable to the D-region or the F-region of 
the ionosphere. The quasi-longitudinal approximation is 
derived from this theory and the range of validity of this 
approximation is illustrated. The restrictions imposed by 
the use of a sharply bounded model ionosphere are discussed. 


Focusing, defocusing, and refraction in a circularly 
stratified atmosphere, K. Toman, J. Research NBS 
4D, No. 3, 289 (1960). 


Focusing, defocusing, astronomical refraction and path length 
of rays as a function of the departure angle A of the ray at 
the source is described for cases with the source outside, 
inside, or on the boundary of a circular stratification. Rel- 
ative to zero elevation angle symmetrical and centrosym- 
metrical distributions are found. 


Impedance characteristics of a uniform current loop 
having a spherical core, S. Adachi, J. Research 
NBS 64D, No. 3. 297 (1960). 


The radiation impedance is derived by the EMF method in a 
convenient form as the sum of the self-radiation impedance 
of a loop in the free space and an additional term due to the 
reaction between the loop and the sphere which is propor- 
tional to the well-known expansion coefficient of a magnetic- 
type scattered wave from a sphere in an incident plane wave. 
rhe first antiresonance frequency has been given in the form 
of a universal curve for a very small uniform current loop 





with core of an arbitrary composition of uw, and e¢,, subject 
to the condition that the refraction coefficient N=yy,«, is 
extremely large. Some numerical calculations show that 
high-z core is desirable for a comparatively lower frequency 
region, and high-e core is rather desirable in an antiresonance 
region. 


Basic theorems in matrix theory, \l. Marcus, NBS 
Applied Math. Series 57 (1960) 15 cents. 


This is a survey of the basic identities and inequalities of 
matrix theory. Included are results dealing with elementary 
properties, canonical forms, invariance, congruence, com- 
mutativity, orthogonalization, eigenvalues, determinants, 
submatrices, rank, determinant and rank inequalities, nu- 
merical methods for inversion and eigenvalues, condition 
numbers. 


Relaxation processes in multistate systems, K. E. 
Shuler, Phys. of Fluids 2, No. 4, 442 (1959). 


An analysis is made of the relaxation of “‘multistate’’ systems, 
i.e., systems with more than two quantum states or two 
different chemical species, for linearized processes described 
by the “Master Equation” of the theory of transport processes. 
These results are compared with those obtained from the 
analysis of two-state relaxation processes and the concept of 
“relaxation time’’ is discussed in this framework. A dis- 
cussion is presented of the transformation of microscopic to 
macroscopic relaxation equations. The existence of periodic 
and/or aperiodic oscillatory solution of the linear multistate 
relaxation equation is investigated. It is shown that the 
multistate relaxation equations admit of aperiodic oscillatory 
solutions. 


Invariant and complete stress functions for general 
continua, (. Truesdell, Arch. Rational Mech. Anal. 
1, No. 4, 1 (1959). 


This paper presents an organized survey of results, with 
proofs, concerning exact general solution of the underdeter- 
mined equations of motion of a continuous medium. Special 
attention is given to a flat space of n dimensions, to an 
n-dimensional space of constant curvature, and to a 2-dimen- 
sional surface. Two methods are used: (1) the representation 
of a solenoidal vector as the curl of a vector potential, and 
2) a variational formulation, essentially the,converse of the 
principle of virtual work. 


The construction of Hadamard matrices, KE. C. Dade 
and K. Goldberg, Mich. Math. J. 6, 247 (1959). 


It is proved that a Hadamard matrix of order 4n can be 
constructed if there exists a transitive permutation group of 
odd order and degree 4n—-1 whose subgroups leaving one ele- 
ment fixed have three transitivity sets each. 


Bending and stretching of corrugated diaphragms, 
R. F. Dressler, ASME Trans. 81D, No. 4, 641 (1959). 


Solutions of the exact linear elastic shell equations for all 
stresses and displacements in a typical corrugated diaphragm 
are presented for three cases over a 9 to 1 range of thicknesses. 
Results were obtained by numerical integration in an elec- 
tronic digital computer. The effect of thickness variation is 
discussed with respect to both stresses and resultants, and 
peak values needed for design purposes are presented,  Cir- 
cumferential and meridional stresses are found to be equally 
important throughout the thickness range analyzed. Bending 
and membrane are likewise equally important 
throughout the range. Peak values in some cases occur near 
the outer rim. 


stresses 
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Modular forms whose coefficients possess multi- 
plicative properties, Mf. Newman, Ann. Math. 70, 
478 (1959). 


Let q be a prime, r and s nonzero integers, and 9 (r) the 
Dodekind modular form. All modular forms of type 
n’(r)n*(qr) such that the related Dirichlet series possesses 
an Euler product are determined. 


The analysis of latin squares with a certain type of 
row-column interaction, J. \landel, Technometries 1, 
No. 4, 879 (1959). 


A serious limitation in the use of latin 
founding of the main effect of each factor 
of the remaining two factors. In some cases, the interaction 
of rows and columns can be expressed as a multiplicative 
term of assigned factors associated with rows and columns. 
The analysis of such designs is presented in detail. An 
example is discussed. The method is extended to the case 
where the treatment effects are linearly related to a given 
concomitant variable. A brief discussion given of the 
relation of the proposed method with other tests for non- 
additivity. 


squares is the con- 
with the interaction 


Is 


A note on algebras, A. J. Goldman, Am. Math. Mo, 
66, 795 (1959). 


The usual sufficient conditions that both regular representa- 
tions of an algebra A be faithful are shown not to be necessary. 
A necessary and sufficient condition (viz, that no nonzero 
member of A be either a left annihilator or a right annihilator) 
is derived. 


Linear differential equations of the second order with 
a large parameter, F. W. J. Olver, J. Soc. Indust. 
Appl. Math. 7, 306 (1959). 


An investigation is made of the differential equation 


dw . 


73 u*p(Z) 


+q(Z)}w 


for large values of the parameter u. 
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